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Editor’s Preface 


Besides ubiquitous radio communication, radar sensor technology is the most 
important application of high-frequency technology today. In addition to the 
classic application fields of radar technology, aerospace and defense, radar 
sensors have recently become a substantial business in the automotive field. 
Although radars in cars today are all still part of pure driver assistance systems 
(i.e. the responsibility lies exclusively with the driver), work on future systems 
for autonomous driving has long been feverish. In autonomous driving, the 
system must take responsibility for any driving maneuvers, which poses a 
major challenge to the reliability of the sensor technology. In addition to the 
increasing resolution and thus complexity of the individual sensors, the tasks of 
the systems are also becoming much more extensive. A very important aspect 
besides the realization of the systems is therefore their validation. Covering all 
critical traffic scenarios in test drives on public roads would, however, require a 
completely unrealistic number of kilometers driven, since the scenarios that are 
critical for safeguarding only occur in a very small percentage of the kilometers 
driven. In addition, the scenarios are not reproducible. For these reasons, 
additional validation methods are indispensable, such as so-called "hardware 
in the loop" (HiL) systems. An essential part of this system is the radar target 
simulator, which must be able to "play" desired target scenarios to a vehicle 
radar, if possible in real time. Radar target simulators have been used for about 
20 years for functional tests of automotive radars, but only in the last few years 
have systems been developed that can also vary the angle of the targets. There 
are mechanical systems and electronic approaches with different shortcomings. 
This is exactly where the work of Dr.-Ing. Axel Diewald comes in. 


In his dissertation, Mr. Diewald has developed important scientific foundations 
for high-precision radar target simulators. The work is based on the system 
concept of a digital radar target simulator implemented on a field-programmable 
gate array (FPGA) with an additional, optional analog target simulator for the 
close range. In addition to basic considerations, the high-frequency and the 
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digital aspects of the realization are considered in the work. For the concept of 
several spatially distributed electronically controllable reflectors, Mr. Diewald 
has, for the first time, fundamentally analyzed the relationship between the 
spatial distribution of the reflectors and the angular resolution of the test radar. 
This results in clear design specifications for future systems. The author’s 
proposed usage of fractional delay filters in the FPGA results in targets that 
can be "shifted" in distance without rasterization. Together with the novel 
concept of lateral target control, atarget simulator is created that allows arbitrary 
and gridless generation of targets in 4D radar target space (3D position plus 
Doppler). Mr. Diewald has succeeded in successfully demonstrating all his new 
concepts in his demonstrator. 


The work of Dr. Axel Diewald represents an important innovation to the state- 
of-the-art and I am sure that his concepts and methods will draw much attention 
and find their way in many applications worldwide. For Dr. Diewald, with his 
creativity, broad technical and great organizational skills, I wish him further 
much success in his scientific engineering career and economic endeavors. 


Prof. Dr.-Ing. Dr. h.c. Thomas Zwick 
- Institute Director - 
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Zusammenfassung 


Automobil-Radarsensoren spielen eine wichtige Rolle bei der aktuellen Ent- 
wicklung des autonomen Fahrens. Ihre Fähigkeit, Objekte auch unter widrigen 
Witterungsverhältnisse zu erkennen, macht sie unverzichtbar für die Umge- 
bungserfassung in autonomen Fahrzeugen. Da ihre Funktionalität direkt am 
Einsatzort validiert werden muss, ist ein vollständig integriertes Testsystem 
erforderlich. Radarzielsimulatoren (RZS) haben in letzter Zeit sowohl in der 
Forschung als auch im kommerziellen Bereich viel Aufmerksamkeit erregt. 
Sie sind in der Lage, ein zu prüfendes Radargerät zu täuschen, indem sie eine 
künstliche Umgebung erschaffen, die aus virtuellen Radarzielen besteht. Dazu 
empfängt ein RZS das vom zu testenden Radar ausgestrahlte Signal, modifiziert 
es im Hinblick auf die gewünschten Zieleigenschaften und sendet es zurück an 
das Radar. Damit diese virtuelle Umgebung so glaubwürdig und realistisch wie 
möglich ist, müssen die virtuellen Radarziele hinsichtlich ihrer Eigenschaften 
möglichst präzise erzeugt werden. 


Die meisten atuellen RZS simulieren die Winkelposition eines virtuellen Ziels 
durch Multiplexen des künstlichen Radarechos zwischen mehreren lateral 
verteilten Antennen. Dies ermöglicht jedoch nur Ziele an diskreten Winkel- 
positionen zu simulieren und führt zwangsläufig zu nicht vernachlässigbaren 
Diskrepanzen zwischen Soll- und Ist-Winkel. In dieser Arbeit wird eine 
neue Methodik vorgestellt, die es ermöglicht, virtuelle Ziele an beliebiger 
Winkelposition zu simulieren, indem die künstlichen Radarechos benachbarter 
RZS-Antennen überlagert werden. Dies ermöglicht nicht nur die durchge- 
hende Simulation von lateral bewegten Zielen, sondern reduziert auch den 
auftretenden Winkelfehler drastisch. 


Die Entfernung eines virtuellen Ziels wird simuliert, indem das eintreffende 
Radarsignal künstlich verzögert wird. Aufgrund der internen Architektur 
derzeitiger RZS kann diese Verzögerung und damit auch die Zielentfernung 
nur in bestimmten Schrittweiten eingestellt werden. In dieser Arbeit wird 
eine weitere neue Methodik vorgestellt, die eine stufenlose Einstellung der 
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Verzögerung und damit einhergehend die Simulation von virtuellen Zielen 
in beliebiger Entfernung ermöglicht. Damit lassen sich sehr kleine Entfer- 
nungsänderungen und somit sogar die Bewegung virtueller Ziele innerhalb 
des Messzyklus eines Radars simulieren, so dass auch nicht-ideale Ziele mit 
realistischen Eigenschaften wie Entfernungs- und Geschwindigkeitsmigration 
erzeugt werden können. 


iv 


Abstract 


Automotive radar sensors play a major role in the current development of au- 
tonomous driving. Their ability to detect objects even under adverse weather 
conditions makes them indispensable for environment-sensing tasks in au- 
tonomous vehicles. Since their functional operation must be validated in-place, 
a fully integrated test system is required. Radar target simulators (RTSs) have 
recently drawn a lot of attention both in research and in the commercial sector. 
They deceive a radar under test (RuT) by creating an artificial environment 
consisting of virtual radar targets. To achieve this, an RTS system receives the 
radar signal emitted by the RuT, modifies it with regard to the intended target 
characteristics and re-transmits it towards the RuT. In order for this environ- 
ment to be as credible and realistic as possible, the virtual radar targets must be 
generated as accurate as possible with regard to their characteristics. 


In most current RTS systems, the angle of arrival (AoA) of a virtual target is 
simulated by multiplexing the artificial radar echo between multiple laterally 
distributed RTS front ends (FEs). However, this only allows the simulation 
of virtual targets at discrete angular positions and inevitably leads to non- 
negligible discrepancies between set and detected AoA. In this work, a new 
technique is presented that enables virtual targets to be simulated at an arbitrary 
AoA by superimposing the artificial radar echoes of adjacent FEs. It not only 
permits the continuous simulation of lateral moving targets, but also drastically 
reduces the occurring angular error. 


The range of a virtual target is simulated by applying an artificial delay to the 
incoming radar signal before re-transmission. Due to the internal architecture of 
current RTS systems, the delay and thus also the target range can only be set in 
certain increments. In this work, yet another new technique is presented, which 
enables a seamless adjustment of the delay and the simulation of virtual targets 
at arbitrary range. This allows to simulate very small range changes and thus 
even the movement of virtual targets within the measurement cycle of a radar 
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enabling the creation of non-ideal targets that feature realistic characteristics 
such as range and Doppler migration. 
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1 Introduction 


1.1 Motivation 


In recent years, the development of advanced driver assistance systems (ADASs) 
and autonomous driving has reached new levels of sophistication. ADASs use 
environment-sensing technologies such as camera, ultrasound, light detection 
and ranging (lidar) and radio detection and ranging (radar) to monitor the sur- 
roundings of the vehicle, detect possible obstacles nearby and react accordingly. 
In order to ensure proper operation and thus the safety of all traffic participants 
involved, autonomous driving functions and the sensors on which they rely 
must be validated thoroughly. Due to their weather robustness and long-range 
capability, radar sensors play a key role in a large share of such systems, and 
must be tested in particular. 


The validation procedure is divided into multiple successive stages and car- 
ried out in parallel with the development process. During the design and 
concept phase, the radar sensor is verified using so-called Model-in-the-Loop 
(MIL) tests; later in the implementation phase, Software-in-the-Loop (SIL) and 
Hardware-in-the-Loop (HIL) tests are employed [FKJ*19]. Finally, during the 
integration phase, field tests are conducted on the road, where the sensor and 
ADAS function are deployed in the vehicle. These field tests are highly labor- 
intensive and require a great deal of effort, as distances in the order of several 
million kilometers must be covered in order to ensure the proper functioning 
of the system [KP16, MGLW15,Sch17b, KW 16]. Moreover, these tests are not 
repeatable, since real-world traffic scenarios are unique, and thus they must be 
reiterated whenever the system undergoes any design changes. 


For this reason, there have been attempts to find test solutions that com- 
plement road tests and thereby reduce the amount of work required. Vir- 
tual vehicle environment (VVE) simulations emulate a digital clone of the 
car and are able to replicate entire traffic scenarios, including the physical 
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and electromagnetic behavior of the vehicle itself as well as its surround- 
ings [FFFG20, DWCP21, SRGD21]. The propagation, reflection and diffrac- 
tion of the electromagnetic waves emitted by the virtual radar sensor are sim- 
ulated using ray tracing technology [AHHD20, Sei21]. Although this method- 
ology represents an unprecedented level of accuracy in virtualized testing, it is 
nonetheless based on software models that can only to a certain degree mimic 
the exact behavior and interactions of real-world components such as antenna 
array and bumper. Hence, sheer software simulations are not realistic enough 
and cannot be used solely for the validation of safety-relevant driving functions. 


AS a consequence, a new testing paradigm is on the rise: Vehicle-in-the-Loop 
(VIL) testing [WHLS19, GMS*18, WMN*20, AVL22]. It enables thorough 
integrative validation of ADAS functions in a controllable laboratory environ- 
ment. As can be deducted from the name, the car is tested as a whole, for 
which a vehicle under test (VuT) is mounted on a roller dynamometer or power 
train test bench that is connected to a VVE simulation [DKK*21]. The torque, 
steering, pitch and roll forces to which the digital car clone is subjected in the 
simulation are applied to the VuT on the test bench through electric motors, 
and its behavioral response is fed back into the VVE simulation. 


In addition to the mechanical-dynamic components, the various sensors of the 
vehicle, including the radar, must also be stimulated and tested. This is one 
of the reasons why radar target simulator (RTS) systems have recently drawn 
a lot of attention, both in research and the commercial sector, as they are 
able to deceive a radar under test (RUT) by creating an artificial environment 
consisting of virtual radar targets. To accomplish this task, RTS systems 
receive the radar signal emitted by the RuT and send back a modified version 
of it. The modifications correspond to the virtual target characteristics that are 
to be simulated, namely range, velocity, radar cross section (RCS) and angle of 
arrival (AoA). They are realized by applying a delay, Doppler frequency shift 
and attenuation to the return signal and controlling the angular position from 
where it is re-transmitted. The basic idea of integrated radar target simulation 
is illustrated in Fig. 1.1. 


1.2 State of the Art 


Virtual Radar Targets < > 


Radar Target 
Simulator 


O» 


Vehicle under Test 


Figure 1.1: Operating principle of integrated RTS-VIL testing 


1.2 State of the Art 


As a matter of fact, RTSs have been around for some time, but have mainly 
been used in other applications, especially in end-of-line testing as part of 
the quality assurance process of radar sensor manufacturers [Roh21a, Key17, 
dSP20, Kon17b, Ros16]. In this context, RTS systems support the calibration 
and functional verification of the sensor after its final assembly and are only 
required to create a single digit number of virtual targets. As the focus now 
shifts towards ADAS validation and VIL testing, the demands on RTSs increase 
drastically as they must now be able to simulate complex traffic scenarios that 
comprise a multitude of radar targets. 


There have been numerous approaches to design multi-target RTS systems, 
using different methods to create the virtual radar targets. The artificial 
radar echo can be generated either by applying a true-time delay [GGS*17, 
EPB16, LEWW14] or a frequency shift [SSGW21, RH20, RISW22] to the 
incoming transmit signal from the RuT, or analyzing & re-synthesizing it 
[DMW18, WMN*20]. In addition, the realization of the selected method can 
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be done using either an analog [EPB16, GGS*17,SJM*21,IMSW19] or digital 
implementation topology [SN20, WMN*20,SP21,JLH21]. It turns out that nei- 
ther implementation concept, analog or digital, is clearly superior to the other, 
as they both have their advantages and disadvantages. The same applies to the 
target generation method, whose selection in combination with the implemen- 
tation topology has major implications on the capabilities and performance of 
the RTS system. 


Creating virtual radar targets with the frequency shift or analysis & re-synthesis 
method is only applicable to radar sensors that employ a chirp sequence modu- 
lation scheme and it requires prior knowledge of the exact waveform parameters 
of the transmit signal. This presents an issue because firstly, this kind of infor- 
mation can only be obtained by extraction, as it is not made publicly available by 
the manufacturers, and secondly, because future radars sensors will potentially 
employ all-digital modulation schemes [GdONA*22, WS14, SGR*21]. There- 
fore, it can be concluded that only the true-time delay target generation method 
meets the requirements of future-proof RTS systems capable of validating au- 
tonomous driving functions. 


Analog RTS systems that utilize the true-time delay method are able to create 
virtual targets in close vicinity to the RuT, but are restricted in their scalability, 
since the generation of each virtual target is assigned to its individual signal 
path within the RTS, which in turn is implemented using dedicated hardware 
components. Evidently, this circumstance prevents any reconfiguration of the 
system at run-time or even post-production and permits the generation of only 
a predefined and limited number of targets. In contrast, digital RTS systems 
are highly flexible due to the reprogrammability of their target generation unit 
and can create a large number of virtual targets. Nevertheless, due to inevitable 
signal delays caused by the conversion to and from the digital domain, the 
minimum simulatable range of a virtual target generated by a digital RTS system 
cannot fall below a certain limit, usually of the order of 20 to 30 m. In order to 
combine the benefits of both implementation topologies, a hybrid system can 
be proposed that splits the incoming radar signal and feeds it into an analog and 
digital RTS in parallel [GSG*18]. Subsequent to the signal modifications of the 
respective target generation units, the signal can be recombined and transmitted 
back to the RuT. More detailed deliberations concerning the implementation 
topologies and the target generation methods can be found in chapter 3. 


1.3 Objectives and Structure of the Work 


The simulation of the angular position of the virtual target is performed inde- 
pendently of the previous remarks. The angle from which the artificial radar 
echo is sent back to the RuT is controlled by multiplexing between multiple 
spatially distributed RTS front ends (FEs) that are arranged in a semi-circular 
formation around the RuT. 


1.3 Objectives and Structure of the Work 


The objective of this work is to design and implement an RTS system that is 
capable of performing sophisticated ADAS validation tests. For this, it must 
not only be able to simulate complex traffic scenarios that consist of a multitude 
of radar targets, but also be compatible with current and upcoming commercial 
radar sensors. As the validation tests are to be carried out integratively and 
in-place, the sensors and their debug interfaces cannot be probed, and no prior 
knowledge of the waveform can be assumed. The system should work with 
all popular modulation schemes and be able to simulate virtual targets that can 
be controlled independently in range, velocity, RCS and AoA. Based on these 
requirements and the previously mentioned considerations, a digital system 
topology and a true-time delay target generation approach were chosen. This 
represents the only solution that not only covers all modulation types, but also 
enables a flexible simulation of a large number of virtual targets. 


As mentioned before, the AoA of virtual targets is simulated by multiplexing 
the artificial radar echo between lateral distributed RTS FEs, which only allows 
for targets at discrete angular positions, inevitably leading to non-negligible 
deviations between set-point and actual AoA. Given the long range capability 
of modern radar sensors, even small angular errors can lead to a preceding 
vehicle being detected in the wrong lane, which can have severe consequences 
for the vehicle occupants. Hence, an RTS system for ADAS validation must 
be designed to minimize the angular error, which is why in this work a new 
technique is presented that enables the generation of virtual radar targets at an 
arbitrary AoA. The basic concept is to superimpose the artificial radar echoes 
of adjacent RTS FEs and control the AoA with the amplitude attenuations of 
the corresponding echoes. This allows to reduce the number of RTS channels 
required and thereby cuts the hardware cost, while at the same time drasti- 
cally decreases the angular error and increases the accuracy with which the 
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virtual targets are created. This new technique overcomes the rasterization and 
discontinuity limitations of the angle simulation when using distributed FEs. 


A similar rasterization applies to the delay that is applied to the incoming 
radar signal. In the digital domain, a true-time delay is induced by buffering the 
incoming samples for a specified number of clock cycles, before they are passed 
on. Since only an integer number of samples can be buffered, the resulting signal 
delay and thus also the target range can only be adjusted in certain increments. 
This determines the range accuracy of the RTS, which must satisfy not only 
the range resolution capability of modern radar sensors, that are on the verge to 
be improved using greater bandwidths, but also the even finer range accuracy 
of the sensors. In this work, yet another new technique is presented, which 
enables a stepless adjustment of the delay so that virtual targets can be created 
at arbitrary range without the restrictions of a rasterization. The application of 
a gridless signal delay is made possible using all-pass fractional delay filters 
whose group delay can be controlled via the filter coefficients. This technique 
also bears the potential to simulate the movement of virtual targets within the 
measurement cycle of a radar allowing to create non-ideal targets that feature 
realistic characteristics such as range and Doppler migration. 


In the following, first the radar fundamentals are explained that are necessary 
for the comprehension of the subsequent elaborations. Next, an overview of 
the state of the art of radar target simulation is given, before the RTS system 
that was built in the course of this work is presented. Thereupon, the two newly 
developed techniques are discussed that overcome the rasterization of range and 
angle, thus enabling the simulation of virtual radar targets at arbitrary positions. 
Finally, a conclusion that summarizes the achievements of this work is given. 


2 Radar Principles 


This chapter addresses the theoretical fundamentals of radar technology in- 
cluding the analog and digital signal processing, which are necessary for the 
comprehension of the following chapters. First, the universal radar range equa- 
tion is derived before subsequently the generic signal model of a chirp sequence 
radar is developed. Next, the digital signal processing steps which enable the 
estimation of the range, velocity and angle of arrival (AoA) of a radar target 
are elaborated. In the end, analytical remarks concerning the resolution and 
unambiguous measurement region of a radar are given. 


2.1 Radar Range Equation 


In order to detect an object with a monostatic radar system, a signal with the 
power of Px is radiated from a transmitter via an antenna with the antenna gain 
G. At the distance R from the antenna the power density is therefore 


PxG 


Stx = (Ar) R2 (2.1) 


The object located at the distance R from the radar sensor reflects the impinging 
electromagnetic wave according to its radar cross section (RCS) a [Sko90]. 
The signal travels back to the sensor, where it is received by the radar with its 
effective antenna aperture of 


Gr? 


Aus’ 2.2 
w= (2.2) 
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resulting in the receive power at the receiver of the radar of [G611] 


_ PuG?oA? 


P= a (2.3) 


The equation is only applicable if the far-field condition is met [Lud13]. 


2.2 Chirp Sequence Radar 


Although the RTS system implemented as well as the methods presented in the 
course of this work are not limited to a specific radar modulation scheme, the 
subsequent analytical descriptions in this work are based on a chirp sequence 
radar since it represents the current de facto standard in automotive radar 
applications [Rob22, Con21]. 


The radar transmits a continuous signal, consisting of multiple chirps, each of 
which has a linear frequency slope. The signal propagates through free space, 
is reflected by a potentially dynamic target and returns to the sensor with a time 
shift r. Thereupon, it is mixed with the transmit signal to form the so-called 
beat signal, which is subsequently low-pass filtered and discretized. The range 
of the target can be estimated by determining the frequency of the beat signal 


Figure 2.1: Chirp sequence radar transmit (blue) and receive (red) signal 


2.2 Chirp Sequence Radar 


which corresponds to the difference of the instantaneous frequencies between 
the transmit and receive signal and hence to the time delay. The potential 
velocity of the target can be estimated by analyzing the relative phase shift 
among the chirps. More detailed explanations concerning the radar signal 
processing can be found in 2.3. An illustration of the transmit and receive 
signal of a chirp sequence radar can be seen in Fig. 2.1. 


2.2.1 Signal Model 


In the following, the signal model of a chirp sequence radar sensor is developed. 
The model is kept generic and will be extended in the subsequent chapter for 
the respective applications. For the sake of simplicity and clarity, it focuses 
primarily on the signal phase in order to facilitate the comprehension of the 
approaches developed in chapter 5 and 6 and their limitations. 


The radar sensor transmits a chirp signal whose frequency increases linearly 
over time and can be described by 


Ri jes z 1, 1€ [0,7] (2.4) 


where fer denotes the start frequency of the chirp, B the bandwidth and T; the 
duration of the chirp. The phase of the signal can by derived by forming the 
integral over the signal frequency 


i B 
xlt) = an f fN- dt =2 | fa tt — -t (2.5) 
0 2T; 
which ultimately enables the deduction of the time signal 
: i B 4 
Stx (t) = exp {jp (t) } = exp 4327 | for- t+ ar t (2.6) 
S 
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neglecting the amplitude. For a single measurement, the radar transmits N. 
chirps, whose phase and corresponding signal can be expressed as 


+ P(r) , ne € [0, Ne - 1] (2.7) 


B 
Yix(t, Wc) = 20 fe -t+ IT. Ta 


Stx (t, nc) = exp {jix (t, Mc) } = exp {i2n B -t+ zZ Ple on) (2.8) 


2T; 


D(n¿) describes the random start phase of the chirp which is caused by resetting 
the frequency back to fer in between the chirps. As stated before, the receive 
signal and its phase can be described as a time shifted version of the transmit 
counterparts 


Six (t, Ne) = Six (t — T, Nc) (2.9) 
Prx(t, Nc) = P(t = T, Ne) (2.10) 


The mixing of the transmit and receive signal, i.e. the de-chirping, can be 
expressed either as the multiplication of the transmit signal with the complex 
conjugate of the receive signal 


sp(f) = S(t, nc) ` Sea: (E Ne) (2.11) 


= exp (on fer + = (27 +t- e)l} 


or as the difference of their respective phases 
polt) = pa (t, ne) — rx (t, Nc) (2.12) 


B 
=2n [fer + (r-e) 


The instantaneous range of a target which is initially located at a distance of 
Ro and features a radial velocity of v throughout a radar measurement can be 
expressed as 


R= Ro + v(t +ncTe) (2.13) 
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where T, denotes the chirp repetition period. Thus, the resulting instantaneous 
delay of the receive signal can be described by 


2R 2v(t + ncTe) 
T= — = 1 + — 
co co 


(2.14) 


where To represents the initial signal delay. Substituting the delay in (2.12) 
yields the beat signal phase of a dynamic target of 


Range-Doppler coupling 


co 
[R(t 8) mE) 16 (E (efo) an 


U 
ee” 
Velocity estimation 


Polt, nc) = 27 


Range estimation 


B2 ) B 2y? 
+22 (1-2}}-2- BP) aay 
Ts co 


(2.15) 


B2v 2v B 
+ (22 | = 22)) -t+ ncTe + To (fr - Žo) 
m c ama 
Range-Doppler migration Angle of arrival estimation 


The first term has a linearly increasing phase along time t, whose slope mainly 
depends on the initial time delay to of the associated target. This can be used 
in the next step to estimate the range of the object. The term also incorporates 
a share that is proportional to the velocity of the target. This part is known as 
range-Doppler coupling and leads to a small error in the range estimation caused 
by a dynamic object, which can be compensated after the Doppler estimation. 
The second term features a linear phase slope along chirp index ne, which is 
related to the Doppler shift 


h= Z fa (2.16) 
co 

that the signal experiences and which can be utilized for the velocity estimation. 
The fifth term takes into account that not only the range of a moving target 
changes during a radar measurement, but also the induced Doppler shift is 
dependent on the instantaneous carrier frequency which in turn changes during 
the transmission of a chirp. This so-called range-Doppler migration plays an 
important role in radar signal processing as will be further elaborated in chapter 
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5. The last term describes the phase offset of the beat signal which will later be 
derived as antenna-element-specific and which can be used for the estimation 
of the AoA of the target when multiple antennas are employed. 


2.3 Signal Processing 


For the subsequent radar signal processing, the beat signal phase from (2.15) is 
simplified with regard to the terms needed for the range and Doppler estimation, 
since only the basic principle of the detection process shall be shown. More 
detailed and mathematically complete derivations can be found in chapter 5.2 
and 6.2. In addition, by taking advantage of the relations 


V X co (2:17) 
B< fo (2.18) 
7) <T (2.19) 


the phase of the beat signal can be further streamlined to 
B 2v 
Yo(t,n.) = 27 T? ‘t+ — fer Nele (2.20) 
s co 


The analog-to-digital converter (ADC) of the radar sensor discretizes the 
beat signal with the sample frequency of f, and thereby substitutes the time- 
continuous variable t with the time-discrete sample index ns 


piers (2.21) 


S 
resulting in a time-discrete beat signal phase of 
B 2 
Yp(Ng, Nc) = 2m | —To Ng + Ay "Ne (2.22) 
Ns A 


where N, = T; fs describes the number of samples per chirp. 
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2.3.1 Range Estimation 


For the range estimation of the target, first, a fast Fourier transform (FFT) is 
applied to the sampled beat signal which corresponds to sample-wise multi- 
plication with a complex exponential term and summation over all samples 
Ns 


Ns-1 
x _. ken 
sr(k, nc) = 2 Sb (Ms, Ne) + EXP (o N. | (2.23) 
2v pa, n 
=exp {ian Fr . nal} . 2 exp [zzz [Bto — k] - a 


where k € [0, Ns — 1] describes the bin index of the newly formed value- 
discrete range spectrum. The expression can be reformulated using the partial 
sum of a geometrical series [BSMMO1] and Euler’s formula [K6n13] to 


Ns-1 à 

s 1 — exp {j27 [Bro — k 

> exp | j2n [Bro k]- "3 | 2 explant Pig Al) (2.24) 
n,=0 Ns 1 —exp {iam [Bto — k] x} 


_ _&xp {j7 [Bro — k]} 
exp (io [Bro — k] x} 
___ exp {ja [Bro = k]} = exp (-j7 [Bro — k]) 
exp {in [Bro — k] a} — exp (ir [Bro — k] a} 


= exp {in [Bm -u (1- 2) See a ED 
sin (x [Bro — k] x) 


S 


For Ir [Bro — k] a < 1, the sine function in the denominator can be substi- 
tuted with its argument, forming a sinc-function 


sin (2 [Bro — k] ) 


sin (x [Bto — k] x) 


= N,- sinc (Bto - k) (2.25) 
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Replacing the sum in (2.23) with these equations allows to simplify the expres- 
sion to 


SR(k, ne) = Ns - sinc (Bro — k) (2.26) 


2v Bto — k 1 
qie mt 1-)|' 


The maximum of the sinc-function can be detected as 


- exp fiar 


kaet = Bro (2.27) 


which in turn enables the estimation of the range of the target to 


C 
Réet = kaet * >. = Ro (2.28) 


2.3.2 Velocity Estimation 


The estimation of the velocity of the target proceeds similar to the range es- 
timation. A second FFT is applied to the output of the first one, resulting 
in 


N.-1 
sp(k,m) = X, sr(k, ne) exp [izr e) (2.29) 


n¿=0 


Bro — k 1 
= N; - exp {i2n | > I -= ral - sinc (Bro — k) 


= 2v n 
. > exp  j2n ram] . a 


n¿=0 


where Tm = T¿N, denotes the complete measurement time and m € [0, N¿— 1] 
describes the bin index of the Doppler spectrum. In the same manner as in 
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(2.23), the expression can be reformulated using the partial sum of a geometric 
series and Euler’s formula to 


Ne-1 I In ET, - 
> ew [720] = Tn -n x) = exp li Haa l (2.30) 
Fm Ne 1 — exp {jr [Tm — m] | 


7 exp {ja [#7 — m|} 
exp {ir [Ta -m]} ~ exp {ie [u = ml} 
“exp [jn [287m = m] a} = exp (im [mu - m] oe] 


= exp fin [Zra -m| (1-2). sin (x [227 |) 


Ne sin (7 [Tm - m] xx) 


Similar to before, the sine function in the denominator can be approximated 


with its argument for Ir [Tm - m] ral < 1, forming a sinc-function 
A 2v 
sin (x [XTn —m]) 2y 
— C z Ne - sine q Im -m (2.31) 
sin (7 [Tm — m] x) 
which together with (2.30) can substitute for the sum in (2.29) leading to 
; . 2v 
sp(k,m) = N,N. - sinc (Bro — k) - sinc Im -m (2.32) 


(Bro - k) I E 7.) (Am m] | 7 x) 


The maximum of the second sinc-function can be detected as 


- exp {in 


2 
Mdet = Im (2.33) 
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from which the Doppler frequency and the velocity of the target can be derived 
to 


Mdet 2v 
=== 2.34 
SD, det Ta 7 (2.34) 
A 
Vdet = Mdet * IT, = Da =y (2.35) 


2.4 Angle of Arrival Estimation 


In order to determine the exact position of the target, in addition to the range 
information, the angle of incidence of the back-reflected electromagnetic wave 
must also be measured. For the case of a single transmit and multiple receive 
antennas, the radar sensor takes advantage of the fact that the impinging wave 
is received at the receive antenna elements with a relative phase offset of 


dix sin(8 
Arn = 2r sin) ies melNe=1 (2.36) 


where d,x describes the lateral distance between the antenna elements, O € 
[-90°, 90°] the incident angle of the reflected wave, n,, the receive antenna 
element index and N,x the number of receive antennas. The phase offset, in 
turn, translates to an element-specific signal delay of 


{ dyx sin(@) 
co 


Nex (2.37) 


T=T0 


which can later be exploited to estimate the AoA. 


2.4.1 Virtual Antenna Array 


In the case of multiple transmit and receive antennas, a virtual antenna array can 
be formed. The principle of a virtual antenna array is based on the simulation of 
a virtual receive array in the digital domain, whose number of channels is greater 
than the number of physically existing transmit (N;x) and receive antennas (N;x). 
The combination of a transmit and a receive antenna forms a channel, the 
corresponding relative phase offset of which behaves as if a larger array of only 
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one transmit but Nix : N;x receive antennas existed. A prerequisite for this is 
that the signals of the different transmitters can be separated unambiguously at 
each receiver. For this, the multiple-input multiple-output (MIMO) method can 
be employed which uses orthogonal or uncorrelated signals for the separation. 
It allows to clearly define a transmit-receive antenna pair to be assigned to each 
of the orthogonal signals in the receiver after their separation, which in turn 
corresponds to one of the elements in the virtual antenna array. This creates a 
virtual antenna array with Na = Nix : N;x antenna elements that is larger in its 
lateral dimensions when compared to the actual, physical array, hence enables a 
better angular resolution and a larger unambiguous angle region [Nus21]. The 
relative phase among the virtual antenna elements is 


da sin(@) 


2 AA» na € [0, Na - 1] (2.38) 


Aga =2n 


where da describes the lateral element spacing and na = nıx + Nex : Mix the 
element index. The phase progression corresponds to an element-specific 
signal delay of 


da sin(@ 
„amd, 
co 


(2.39) 


T=T0 


2.4.2 Digital Beamforming 


In order to extract the angle information from the relative phase offset among 
antenna elements of the impinging radar wave, digital beamforming can be 
applied. This process operates similar to a Fourier transform insofar that 
it creates a (pseudo) spectrum that enables to determine the signal’s spatial 
frequency. In this case, it is a spatial spectrum that results from the lateral 
distribution of the antenna elements and which correlates with the incidence 
angle. Digital beamforming can theoretically be applied either before or after 
the range-Doppler processing, as it is compatible with time and frequency 
sampling, but is often performed after the range-Doppler detection to save 
precious hardware resources [Sit17]. 


For the subsequent processing steps, the phase of the beat signal from (2.15) is 
adjusted to take into account the antenna-element-specific delay in (2.39). As 
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only the basic principle of the digital beamforming process shall be shown, the 
expression is reduced to the angle-dependent terms 


enor | B da sin(0) a B E sin(@) m) 


Jes Is co 2T; co 


B 
+70 fa = =) (2.40) 


which in turn can be further simplified using the relations in (2.18) and (2.19) 
to 


Yp(na) = 20 u “nat mofa] (2.41) 
sp(na) = exp {jyn(na)} (2.42) 


To estimate the AoA of the target, the one-dimensional beamsteering vector 


ne 
faso) na) (2.43) 


bıp(na,a) = exp (2 
where a € [-90°, 90°] describes all evaluation angles to be considered and 
can be discretized arbitrarily. The vector can be used to apply a Fourier 
beamforming to the beat signal from (2.42) 


Na-l 
2 d r 

sala) = Y) s(a) exp | jn RO na) (2.44) 

__Nacl 

MAT 
a  dasin(a) M-1 
= > sp(na) < exp 4 —j27 na = 
na=0 A 2 


da sin(@) Na - |) 


exp fiar [ofe + 1 > 


Na-1 
. 2, exp {i2n 


na=0 


(sin(9) — sin(a)) -na 


| 
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The remaining sum can, as in the cases of the range and velocity estimation, be 
reformulated using the partial sum of a geometric series and Euler’s formula to 


Na-1 
y exp fiar 


na=0 


| (2.45) 


a (sin(@) - sin(a)) - na 


1 — exp ar [sin(@) — sin(a)] Na) 


1 — exp {jana [sin(0) — sin(a)]) 
m {ins [sin(@) — sin(a)] Na) 


exp lin [sin(@) — sin(a)]) 


na [sin(@) - sin(a)] Na) — exp (rd [sin(@) — sin(a)] Na) 
exp {ix nú [sin(@) — sin(a)]) — exp (inde [sin(@) — sin(a)]) 
sin (x na [sin(@) — sin(a)] Na) 


sin (z4 [sin(@) — sin(a)]) 


exp {in 


= exp {in [sin(@) — sin(a)] (Na - »| . 


As before, the sine function in the denominator can be substituted with its 
argument for rá [sin(@) — sin(a)]| < 1, forming a sinc-function 


sin (x n% [sin(8) - sin(a)] Na) 


sin (x na [sin(6) — sin(a)]) 


= Na sinc (> [sin(@) — sin(a)] ns] 


(2.46) 
which together with (2.45) can replace the sum in (2.44) resulting in 
. (da,. f 
sala) = Na: sinc > [sin(@) — sin(a)] Na (2.47) 
da sin(@ -1 
- exp {ian rf u ) a | 
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From this, the AoA of the target can be detected as 


sin (get) = sin (0) (2.48) 
Uder = 0 (2.49) 


Two-dimensional Digital Beamforming 


Modern radar sensors tend to have a two-dimensional distribution of their 
transmit and receive antennas, which enables them to resolve targets in the 
azimuth and elevation domain simultaneously. The horizontal and vertical 
positions of the virtual antenna elements can be described exemplarily as follows 


Paz = daz,tx “Nyx + daz,rx * Arx (2.50) 
Pa= del,tx “Nx + del,rx “rx (2.51) 


where d,z/el,tx/rx denotes the horizontal/vertical spacing of the transmit/receive 
antennas, respectively. With this, the relative phase offset of the virtual receive 
antenna elements can be described as 


Aga =2 7 


(2.52) 


where y € [-90°,90°] is the incidence angle of the reflected wave in the 
elevation plane. This in turn results in an element-specific signal delay of 


sin(@) cos(W) - Paz + Sin(W) - Pel 
To + co 


Te 


(2.53) 


For the subsequent processing steps, the phase of the beat signal from (2.15) 
is adjusted to take into account the two-dimensional, antenna-element-specific 
delay in (2.53). Similar to before, only the basic principle of the digital beam- 
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forming process shall be shown, which is why the expression is reduced to the 
angle-dependent terms 


Yu(na) = 27 (2.54) 


B 
Jfa- T 


S 


) sin(@) cos(y) "Paz + sin(y) * Pel 
co 


B_(sin(9) cos() + Paz + sin() + pa | B 
HT. (PO pu tae +70 (fs = = 


Considering that in conventional radar sensors the transmit and receive antennas 
are each distributed in only one spatial dimension [Sys21, Har14], the position 
of the virtual elements can be simplified to 


Paz = tx * Mx (2.55) 
Pel = Ax ` Mx (2.56) 


and thus the beat signal and its phase, also taking into account the relations in 
(2.18) and (2.19), can be reduced to 


OO eee | 


Sp (Mtx, yx) = exp {jp (na) } (2.58) 


Pb (Nx, Mrx) = 27 


With the two-dimensional beamsteering vector 


bop (Nix; Mrx, Y, B) = exp [2 [dx sin(@) cos(ß) + mx + dix sin(ß) - ns 
(2.59) 


21 


2 Radar Principles 


a two-dimensional Fourier beamforming can be applied, leading to 


Nix-1 Nrx-1 
2 2 


sala D= Y È, oC m) > bap (Mx, Mx, 0, B) (2.60) 


Nix-1 Nrx-1 
n=- 5 n=- 


x=1 d -1 
ii [er poses © sin(a) cos(8) + sin), | 


N-i Nx-1 27 
y exp (i= |an (sin(@) cos(yw) — sin(a) cos(8)) + nrx 


Nx=0 n,x=0 
+ dix (sin(w) — sin(8)) - n|} 


The two sums can be simplified individually utilizing the same procedure as 
before, resulting in 


sala, B) = Na : exp {ian rf + = sin(@) cos(W) M! (2.61) 


- sinc [E [sin(@) cos(yw) — sin(a) cos(ß)] No] 


d Nx-1 
e sin(w) = 


- sinc (= [sin(yw) — sin(ß)] Na) 


From this, the azimuth and elevation AoA of the target can be detected as 


sin (Ager) cos (Baer) = sin (0) cos (y) (2.62) 
sin (Baer) = sin (y) (2.63) 

Qdet = 0 (2.64) 

Bact = Y (2.65) 


22 


2.5 Limitations 


2.5 Limitations 


Depending on the parameterization of a radar sensor, different limitations arise 
with regard to its achievable resolutions and unambiguous measurement re- 
gions. 


2.5.1 Resolution 


In radar applications, it is important to distinguish between resolution and 
accuracy. The resolution of a radar describes its ability to discriminate two 
objects and cannot be improved through sheer signal processing, whereas the 
accuracy represents the precision with which the parameters of a target can be 
estimated which can be improved, for example, by zero padding or interpolation 
methods [RSSH 10]. 


Range Resolution 


The range resolution of a radar specifies how close together two point targets 
with identical velocity and RCS can be in order to still be perceived as two 
separate objects. It follows from the consideration that the resolution of the 
beat signal frequency is equal the inverse of the sampling time [PL15] 


1 
Ah= = (2.66) 


which together with the relation of the constant frequency slope of 


B_f 
=== 2.67 
e (2.67) 
yields the range resolution of 
co 
AR = — 2.68 
>B (2.68) 
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Velocity Resolution 


Similar to the range resolution, the definition of the velocity resolution follows 
the consideration that the resolution of the Doppler frequency is equal to the 
inverse of the measurement time 


1 
Ab = z- (2.69) 


which can be substituted in (2.16) to give the velocity resolution of 


co 


Av = >>> 
2 ferTm 


(2.70) 


Angular Resolution 


Proceeding in the same manner as before, the spatial frequency along the spatial 
sample points, i.e. the antenna elements, can be set in relation to the AoA of the 
target by deriving the relative phase in (2.38) according to the spatial sample 
index 

1 OPA sin(@) 


fa = In dan) a (2.71) 


The frequency resolution is again defined by the inverse of the (spatial) sampling 
period 


1 


Afa = 2.72 
fa EA (2.72) 
This allows to express the analytical angular resolution as 
AQana = arcsin | (2.73) 
ANA 


All of the above definitions for the respective resolutions are to be regarded as 
theoretical statements only. Due to hardware imperfections and the influence of 
noise, they cannot be achieved in practice. For the range and velocity resolution, 
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this circumstance is not decisive, since they are quite fine and factorial deviations 
have only a negligible impact. 


Of greater importance, however, is the precise estimation of the angular resolu- 
tion, which is why two adapted definitions can be found in literature. Both are 
based on considerations concerning the sinc-function that arises after the beam- 
forming process in (2.47). The first one states that two targets can be resolved 
in the angular domain, if the maxima of the corresponding sinc-functions are 
separated by at least twice the half-power beamwidth [Bal15]. The half-power 
beamwidth, in turn, can be found by solving 


sinc (Xhp) = > (2.74) 


numerically which results in xpp = 0.443. With this and sin(x) ~ x, the first 
literature angular resolution definition can be expressed as 


0.8864 


_—— 2.75 
TAN, (2.75) 


AQhp = 
The second definition is based on the Rayleigh criterion and determines the 
angular resolution for the case that the maximum of the first sinc-function 
coincides with the first minimum of the second [Hec12, HTS*12] which leads 
to 


Nye = ——— (2.76) 


A common constraint for both definitions is that the two targets to be resolved 
are incoherent in regards of the phase of their corresponding radar echoes. 
This means that the superposition of their reflected radar signals occurs only 
in accordance to the signal powers, not the phases [RSSH10]. Moreover, the 
definitions are only accurate for targets at the antenna array’s broadside at 9 = 0° 
and worsen with increasing AoA, since the sinc-function broadens. 
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2.5.2 Unambiguity 


The unambiguity regions indicate up to which maximum range, velocity or 
angle of an object the radar sensor can still estimate its parameters correctly. 
Beyond this limit, the target can still be detected — except for the AoA - if the 
sensor is not limited by noise, but then the estimated value no longer corresponds 
to the actual one. This is ultimately due to the fact that if the sampling theorem 
is violated, the detection which lies above the maximum value is folded into the 
previous unambiguity region. 


Maximum Unambiguous Range 


The maximum unambiguous range is defined by the maximum frequency of 
the beat signal that can still be measured unambiguously, which in turn is given 
by the Nyquist-Shannon sampling theorem [Whi15] 


fona = > (2.77) 


(2.78) 


Maximum Unambiguous Velocity 


According to the sampling theorem, the maximum unambiguously measurable 
Doppler frequency is defined by the chirp repetition period 


1 
JD,ua = oT. (2.79) 
which leads to the maximum unambiguous velocity of 


co 


= —— 2: 
Vua AT fos (2.80) 
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Maximum Unambiguous Angle of Arrival 


Applying the sampling theorem to spatial sampling, the maximum unambigu- 
ously measurable spatial frequency (2.71) is defined by the spatial sampling 
period da [Har14] 


1 


FA,ua = Od, (2.81) 


which results in the maximum unambiguous AoA of 


(2.82) 


a 
Qya = arcsin (5) 
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This chapter elaborates on the state of the art of radar target simulation, both 
in research and the commercial sector. First, the basic principle as well as 
the different modes of operation of RTS systems are discussed, before their 
miscellaneous applications and use cases are presented. Next, the distinctive 
characteristics of a radar target and the corresponding signal modifications 
that an RTS system must perform in order to simulate virtual radar targets are 
explained. The two possible implementation topologies to realize these signal 
modifications are presented before finally the topologies as well as the different 
concepts for the target generation and the angle simulation are compared. 


3.1 Basic Principle 


The basic idea of an RTS system is to deceive a radar under test (RuT) by 
creating an artificial environment that comprises of one or multiple virtual radar 
targets. For this purpose, the RTS generates synthetic radar signal echoes and 
emits them towards the RuT which interprets them as reflections of real-world 
objects. In order to create these artificial echoes, the RTS system receives the 
radar signal emitted by the RuT, down-converts and modifies it with regard to 
the desired target characteristics and subsequently up-converts and re-transmits 
it towards the RuT. The modifications of the radar signal can be realized using 
an analog or digital system topology as is further explained in chapter 3.4. 
There are alternative solutions where the incoming signal is not modified, but 
analyzed and re-synthesized with alterations that correspond to the desired 
target properties. The operating principle of both concepts are discussed in 
more detail in chapter 3.3 and their pros and cons are elaborated in section 
3.9.2; 
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3.2 Applications 


Radar target simulation is mainly used in three different applications. The first 
one is end-of-line production testing as part of the quality assurance process of 
radar sensor manufacturers. It serves the purpose of calibrating the sensor and 
verifying its functionality after its final assembly. The RTS systems developed 
for this application [Roh21a, Key17, dSP20, Kon17b, Ros16] simulate only a 
single target at high accuracy and the processed output of the radar sensors is 
set in reference to the ground truth. 


The next use case is the development of ADAS functions that are implemented 
directly on the hardware of the radar sensor. In this context, RTS systems 
provide a powerful debugging tool that enables the stimulation of the RuT 
with synthetic radar echoes while monitoring the output of the driving function 
under development [Key22, dSP21]. In this application, the sensor is tested 
independently and segregated from its later application site. 


In contrast, the third use case concerns the in-place verification of the radar 
sensor integrated in an autonomous vehicle including the interaction with its 
connected devices. For this, a vehicle under test (VuT) is placed on a test bench 
that is wired to a virtual vehicle environment (VVE) simulation in which a 
digital clone of the car is emulated. The torque, steering, pitch and roll forces 
that the clone is exposed to are applied to the VuT on the test bench through 
electric motors and its behavioral response is fed back into the VVE simulation. 
The integrated radar sensor of the VuT is stimulated by an RTS system which 
receives a list of the virtual radar targets that it is supposed to simulate from the 
VVE simulation. RTS systems that are employed for this application must be 
able to generate a multitude of different synthetic radar echoes simultaneously 
in order to simulate complex traffic scenarios and thus suffice the requirements 
of ADAS validation tests [Roh21b, WMN*20, GMS* 18, BAB*21, GSG*18]. 
This setup is considered a Vehicle-in-the-Loop test system and enables holistic 
verification of autonomous driving functions with the radar sensor deployed in 
its actual operation site. The basic concept of the setup is illustrated in Fig. 3.1. 
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Figure 3.1: Basic concept of an RTS system integrated in a Vehicle-in-the-Loop test bench for the 
application of ADAS validation testing 
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3.3 Target Generation Methods 


As discussed in chapter 2, a radar target is characterized by four distinct prop- 
erties: range, velocity, electromagnetic reflectivity (cf. RCS) and angular po- 
sition. These characteristics correspond to alterations of the received signal of 
the radar sensor in comparison to the transmitted one in terms of delay, Doppler 
frequency shift, amplitude and phase offset among the antenna elements of the 
RuT. Hence, an RTS system must simulate all of these signal amendments in 
order to generate a virtual radar target. 


3.3.1 Range, Doppler and RCS Simulation 


As mentioned before, an RTS system creates virtual targets by receiving the 
emitted radar signal from the RuT and sending back a modified version of it. 
The altered signal can be generated either by modifying the incoming original 
radar signal [IRW21,Roh21b, GGS* 17] — hereafter referred to as the modifica- 
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tion approach — or analyzing and re-synthesizing it with changes [DMW18] — 
hereafter called the re-synthesis approach. 


The Modification Approach 


The modification approach, in turn, can be divided into two different imple- 
mentations, which differ only in the simulation of the range of the target. In the 
first implementation, this target property is mimicked by delaying the incoming 
radar signal by the same amount as it would be if it were propagating in free 
space to and from the real-world object [EPB16,GGS*17, LEWW14]. Thus, 
there required latency that the RTS needs to apply is 


Trts = Am (3.1) 
co 
where Ryim describes the range of the simulated target. The second implementa- 
tion takes advantage of the fact that FMCW and chirp sequence radar receivers 
transform signal delay into a beat frequency. Hence, shifting the frequency of 
the incoming radar signal at the RTS by 


2BRsim 
Rus = g, s = 


2 
Toco (3.2) 


enables to control the frequency of the beat signal directly, which is later 
interpreted by the radar sensor as the target range [RH20, RISW22, SSGW21]. 


The simulation of the RCS and the velocity of the virtual target is the same for 
both implementations. The latter is realized through a frequency shift of 


2 ferVsim 
co 


JD.rts = (3.3) 


which corresponds to the Doppler shift that the signal would experience if it 
were reflected from a dynamic real-world target. 


The target’s RCS is simulated through attenuation or amplification of the in- 
coming radar signal. The signal power at the receiver of an RTS system at a 
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distance of R. from the RuT and with a receive and transmit antenna gain of 
Gre is 
_ PxGGkGrsA? 
Pris,rx = (APR (3.4) 
Within the RTS an artificial power gain G,ı is applied and subsequently the 
signal is re-transmitted towards the RuT where the signal power at the receiver 
can be expressed as 


_ PxG?Ge Gua? 


rer (3.5) 


The required artificial gain to simulate a virtual target with an RCS of Grim at 
a distance of Rsim can be derived when equating (3.5) with the receive power 
caused by a real-world target from (2.3) 


Osim 4TR* 
4 2 
Rsim Ge a? 


Gris = (3.6) 


The Re-Synthesis Approach 


The re-synthesis approach requires to analyze the incoming radar signal and to 
extract certain signal parameters, namely the chirp slope B/T, and repetition 
time T, as well as the start and stop frequency. Based on these, a new chirp 
signal is synthesized which is altered with regard to the frequency and amplitude 
modifications that are needed to simulate the velocity and RCS of a virtual target. 
The created artificial signal can also be a superposition of many virtual targets. 
The signal is synchronized in time and periodicity with the incoming radar 
signal and transmitted towards the RuT with a delay that corresponds to the 
range of the target. 


3.3.2 Angle Simulation 


The simulation of the angular position of the virtual radar target is indepen- 
dent of the range, velocity and RCS simulation approaches discussed above. 
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It can be realized either by mechanically rotating the RTS system centric 
around the RuT [ABB*21, dSP21], by rotating the sensor itself [Kon17a], or 
by electronically switching among multiple, spatially distributed RTS anten- 
nas [GMS* 18, Key22, Roh21b]. The advantages and disadvantages of the three 
different concepts are discussed in section 3.5.3. 


3.4 System Topology 


The signal modifications required to generate virtual radar targets can be im- 
plemented in the analog or digital domain. The selection of the topology of 
the target generation process determines certain features of the respective RTS 
system such as the minimum simulatable target range, the system’s scalability 
and the hardware costs. Furthermore, not every target generation method pre- 
sented in section 3.3.1 can be implemented in either, analog or digital, domain, 
which presents an additional consideration to be made when designing such a 
system. Only the simulation of the AoA of the virtual target is independent of 
the selected system domain. In the following, both topologies as well as their 
possible realizations are presented. In chapter 3.5.1 their pros and cons are 
listed. 


3.4.1 Analog Target Generation 


The analog target generation employs discrete RF components in order to apply 
the required signal modifications. Using this topology, only the aforementioned 
modification approach with either a true-time delay [EPB16,GGS*17] or a 
frequency shift concept [SJM*21, IMSW19] for the simulation of the target 
range can be implemented. Block diagrams illustrating the basic design of an 
analog true-time delay and frequency shift RTS system are shown in Fig. 3.2. 


When implementing the true-time delay concept for the simulation of the target 
range, so-called delay lines can be employed. The idea is to delay the signal by 
passing it through either a set of optical [EPB16, LEWW 14] or coaxial cables 
[GGS*17] whose respective combined length correlates with the intended target 
range. The architecture of the delay line set can be designed in one of two ways. 
Either with segments of different lengths, that correspond to a power-of-two 
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RuT 


RuT 


(b) Frequency shift concept 


Figure 3.2: Block diagram of an analog (a) true-time delay and (b) frequency shift RTS system 


series and which the signal either passes through or bypasses [LEWW 14], or 
with segments of equal length, after each of which the signal can be branched off, 
with a following switch matrix to de-multiplex the multiple outputs [GGS* 17]. 
Alternatively, a true-time delay can be induced with surface acoustic wave 
(SAW) filters [ALRP*17]. For the frequency shift concept, a complex-valued 
mixing scheme can be employed, which is realized either through an in-phase 
quadrature (IQ) mixer or an analog multiplier [SSGW21, RH20, IMSW 19]. 
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The analog simulation of the velocity of a virtual target can be performed 
either by direct frequency mixing using an IQ mixer or an analog multiplier 
[IRW20,GMS*18,AOHH21], similar to the frequency shift range simulation, or 
by introducing a slight frequency offset between the up- and down-conversion 
prior respectively subsequent to the other signal modifications [KHN*21]. 


The simulation ofthe RCS of atarget is accomplished with variable gain ampli- 
fiers (VGA) [GGS* 17] that are placed in the signal path and whose amplification 
is adjusted according to (3.6). 


3.4.2 Digital Target Generation 


For the digital generation of virtual radar targets, the incoming radar signal 
is digitalized prior to the application of the required modifications using an 
analog-to-digital converter (ADC). The digital signal samples are fed into 
a field programmable gate array (FPGA) or digital signal processor (DSP), 
which perform the signal amendments [SN20] or the signal analysis and re- 
synthesis [WMN*20] in accordance to the target characteristics to be simu- 
lated. Subsequently, the signal is converted back to the analog domain using 
a digital-to-analog converter (DAC) and re-transmitted towards the RuT. The 
block diagram of a digital RTS system is shown in Fig. 3.3. 


Figure 3.3: Block diagram of a digital RTS system 
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The digital range simulation of a virtual target is achieved through sample 
buffering, storing the incoming radar samples in memory and reading them 
back with a controllable temporal offset [SP21]. Alternatively, the range can be 
simulated with the re-synthesis approach by synchronizing the newly created 
artificial signal with the incoming radar signal in time and controlling the time 
shift with which it is sent back to the RuT [DMW18]. This requires a trigger 
mechanism that detects the arrival of a chirp to which the re-synthesized signal 
is then aligned. 


The velocity simulation of a target is performed, in the case of the modification 
approach, with a digital complex-valued mixer that is often implemented as a 
numerical multiplier [JLH21]. The Doppler signal, that the incoming radar sig- 
nal is mixed with, is generated using a direct digital synthesis (DDS) component 
that is capable of synthesizing signals of arbitrary frequencies [KBS*22]. 


The RCS of a virtual target is simulated either by a simple digital multiplication 
of the radar signal with an adjustable factor [SN20], or by digitally controllable 
VGAs in the analog signal path following the DAC [MRP07]. The two options 
can be combined which yields an extended dynamic range of the system. 


For the re-synthesis approach, the simulation of the velocity and RCS of a 
virtual target is taken into account during the generation of the artificial radar 
signal [WMN*20]. 


3.5 Trade-offs 


In the following, the advantages and disadvantages of the different system 
topologies, target generation approaches, and angle simulation concepts are 
elaborated. 


3.5.1 System Topology 


The main differences that set analog and digital RTS systems apart are the 
minimum simulatable target range and the scalability. Due to the implemen- 
tation of the target generation process on an FPGA or DSP, digital RTSs are 
highly scalable. Additional virtual targets can easily be simulated through the 
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reconfiguration of the platform. The amount and distribution of target gener- 
ation units is only limited by the available hardware resources on the FPGA 
or DSP employed. However, this topology is associated with the disadvantage 
that the radar signal experiences an inevitable delay caused by the ADC and 
DAC, which in turn leads to an unavoidable target range offset. Even the usage 
of low-latency converters results in a minimum target range in the order of tens 
of meters, which can be severe for the entire setup, depending on the use case. 
In contrast, analog RTS systems do not deploy such converters, hence very 
small target ranges can be achieved. Nonetheless, analog setups are not as scal- 
able since the signal modifications needed for the target generation are realized 
with discrete electrical components, thus the entire analog signal path must be 
replicated in order to insert additional targets, which cannot be performed at 
run-time or even post-production. 


Other differences between the two topologies concern the implementable target 
generation methods, their monetary cost, and the physical space consumption. 
Analog RTSs can only realize the modification approach, that is the true-time 
delay and the frequency shift concept. The analysis and parameter extraction, 
that is needed for the generation of an artificial radar signal in the re-synthesis 
approach, requires signal processing steps that can only be performed in the 
digital domain [DMW18]. With a digital RTS system, on the other hand, all 
of the target generation methods mentioned above can be realized. Another 
consideration to be made is that the acquisition cost of a high-performance 
FPGA or DSP, along with high-speed ADCs and DACs, exceeds the price of 
the relatively cheap electrical components of an analog RTS. However, the latter 
is by far more space consuming due to the large spatial extent of the delay lines. 


In order to combine the benefits of both system topologies, there have been 
efforts to build a hybrid system that splits the incoming radar signal after 
the down-conversion and feeds it into an analog and digital RTS in parallel 
[GSG*18]. Subsequent to the signal modification of the respective target 
generation units, the signal is recombined and transmitted back to the radar. 


3.5.2 Target Generation Methods 


The selection of the target generation methodology, that is the re-synthesis and 
the modification approach, the latter of which in turn can be divided into the 
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true-time delay and the frequency shift concept, has major implications on the 
capabilities and performance of the RTS system. In the following, these three 
target generation concepts are compared with regards to system topologies that 
can be used to realize them, processable radar modulation schemes, scalability, 
and minimum simulatable target range. 


Analysis & Re-Synthesis 


The re-synthesis approach does not have any limitations regarding the minimum 
target range that can be simulated. Once the generated artificial signal is syn- 
chronized in time with the incoming radar signal, its transmission towards the 
RuT can be performed with an arbitrary time shift, which can even be negative. 
Furthermore, through the superposition of multiple generated artificial signals, 
a multitude of virtual radar targets can be created whose number is only limited 
by the digital hardware resources of the implementing system [WMN*20]. On 
the downside, the approach can only be realized on a digital RTS setup, as men- 
tioned before. Additionally, it only complies with FMCW or chirp sequence 
radar sensors, which have a constant pulse repetition period and an invariant 
frequency slope. This might present an exclusion criterion for future-proof RTS 
systems, since upcoming radar sensors will possibly employ all-digital modu- 
lation schemes [GdONA*22, WS14,SGR*21]. Moreover, at least the first chirp 
of each measurement is lost, since it is needed for analysis and the extraction of 
the chirp parameters for the subsequent signal re-synthesis. Lastly, the velocity 
simulation of a target poses a problem, since the phase at the beginning of 
each chirp is random, which cannot be compensated by the generated artificial 
signal. 


Frequency Shift 


The frequency shift approach is also not limited regarding the minimum sim- 
ulatable target range, since the utilized frequency to alter the incoming radar 
signal can be chosen freely. Scalability does not present an issue either, because 
multiple frequency shift alterations can be applied simultaneously through an 
analog multiplier. The number of targets that can be simulated is again only 
limited by the digital hardware resources that create the mixing signal. Even 
though the approach can be implemented using either system topology, only an 
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analog realization makes sense, since the conversion to and from the digital do- 
main is associated with the aforementioned drawbacks and in this case a digital 
implementation offers no advantages over the former. Like the previous target 
generation approach, the frequency shift concept is only applicable to FMCW 
and chirp sequence radars and also requires knowledge of the chirp parameters 
a priori. In addition, the concept is susceptible to the unintended generation of 
ghost targets that result from non-idealities in the mixing process. IQ imbal- 
ances cause a frequency shift in the negated direction and thus create targets 
with negated range offsets. Furthermore, intermodulation products cause ghost 
targets at multiples or combinations of the intended target ranges. 


True-Time Delay 


In contrast to the two aforementioned methods, the true-time delay target gen- 
eration concept is able to operate with radar signals of any modulation scheme. 
The signal modifications performed are in accordance with the alterations that 
the emitted radar signal would undergo in reality and therefore are not specific 
to FMCW or chirp sequence. The concept can be implemented using either 
system topology, but each entails limitations. The analog-to-digital and sub- 
sequent digital-to-analog conversion required for the digital implementation 
introduces an inescapable signal latency and thus determines the minimum 
simulatable target range. The use of low-latency converters can mitigate but 
not completely avoid this effect. In the analog implementation, the scalability 
is constrained, since the signal modifications are performed in distinct signal 
paths that employ discrete electrical components and must be replicated for 
each additional virtual target. Finally, the approach is subject to a rasterization 
of the simulatable target range, which is driven by the smallest delay increment 
of the target generation. The limitations of this particular handicap can be 
overcome in digital RTS systems by a new technique presented in chapter 5. 


3.5.3 Angle Simulation 


In the following, the advantages and disadvantages of three different target angle 
simulation concepts are discussed. For the first concept, the RuT is mounted 
on a gimbal, that is a motorized cardanic suspension, and the inclination of 
the sensor is changed in accordance to the target angle to be simulated. This 
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presents a high-precision and low-cost solution, but is, however, only capable to 
simulate a single independently moving angle. Therefore, the concept is only 
applicable to end-of-line testing of radar sensors, but not to the validation of 
ADAS functions. 


In the second concept, the RTS system, or more specifically, its receive and 
transmit antennas, are moved centrically around the RuT using a mechanical 
sledge traveling in a semicircle. This also allows a high-precision simulation 
of the target angle and furthermore enables the realization of multiple indepen- 
dently laterally moving virtual targets. However, due to the possible mutual 
obstruction and occlusion, the number of targets is limited. Moreover, the travel 
speed of the sledge constraints the lateral velocity of virtual targets. 


For the third concept, the receive and transmit antennas of the RTS system are 
distributed statically in a hemispheric formation around the RuT. The lateral 
movement ofa virtual target is realized through electronically switching among 
the antennas. This represents the only of the three solutions capable of simu- 
lating complex traffic scenarios that consist of a large number of independently 
moving radar targets, as required for the thorough validation of autonomous 
driving functions. A major disadvantage of this concept is that virtual targets 
can only be simulated at the discrete angular positions at which the antennas of 
the RTS system are located. This restriction can be tackled with anew approach 
presented in chapter 6 that expands on the concept of spatially distributed RTS 
antennas and enables the simulation of radar targets at arbitrary angles. 
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In this chapter the digital RTS system that was designed and implemented in 
the course of this work is presented. The system allowed to verify the concepts 
developed in this thesis, analyze their limitations and investigate unexpected 
phenomena. In the following, first, the basic structure of the test system, 
consisting of modular front ends (FEs) and a common back end based on a 
field programmable gate array (FPGA), as well as their individual components 
are explained. Hereafter, a link budget and a spectral model of the system are 
developed, which serve to illustrate and reason the design process. Next, the 
target generation process and the implementation of the signal modifications 
included therein, namely delay, frequency shift, phase shift and attenuation, 
are elaborated. Subsequently, a balance of how many resources of the FPGA 
were used for the implementation of the target generation unit is drawn, and an 
estimate of the scalability of the system with respect to the utilized platform 
is given. At the end of this chapter, a generic signal model is developed that 
analytically describes the propagation of the signal from the RuT, through the 
RTS and back and allows the calculation of the modification parameters required 
to simulate targets with desired properties. 


4.1 Basic Design 


The basic design of the realized RTS system consists of modular, individual 
FEs and a common digital back end and is shown in Fig. 4.1. As can be 
seen, the FEs are separated into dedicated receivers and transmitters, whereas 
the back end combines multiple analog-to-digital (ADC) and digital-to-analog 
converters (DAC), as well as target generation instances, and therefore is able 
to serve multiple FEs. The composition of one receiver front end, ADC, 
target generation instance, DAC, and transmitter front end is considered an 
RTS channel. The modularity of the design and the separation of front and 
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Figure 4.1: Basic design of the digital RTS system 


back end enables an easy extension of the system towards a larger number 
of RTS channels and the free configuration of spatial distribution of the FEs. 
Moreover, the common back end facilitates the synchronous adjustment and 
mutual calibration of the individual channels, as required for the concepts 
presented in the work. 


The receiver FE picks up the radar signal emitted by the RuT and down-converts 
it to an intermediate frequency fits, so that it can be digitalized by the ADC of 
the back end. The signal then passes through the target generation unit on the 
FPGA where the desired target properties are applied. The unit is controlled 
by a software application that is hosted by the processing system (PS) of the 
back end, which receives the list of targets to be generated by either a virtual 
environment simulation or user input. The output of the target generation unit 
is fed into the DAC which synthesizes the modified radar signal and passes it 
on to the transmitter FE, where it is up converted back to its original carrier 
frequency and re-transmitted towards the RuT. Front and back end of the system 
are connected through coaxial cables. 


4.1.1 Front End 


The modular FEs are divided into dedicated receivers and transmitters that 
receive and down-convert the radar signal emitted by the RuT or up-convert 
and re-transmit the signal modified by the RTS. Block diagrams of both re- 
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Figure 4.2: Block diagram of the (a) receiver and (b) transmitter front end 


ceiver and transmitter are shown in Fig. 4.2 with reference indications to the 
spectral model in section 4.3. The LO signal that is used for the down- re- 
spectively up-conversion of the radar signal passes through several frequency 
multiplication stages which allows for the signal to be fed to the FE at a much 
lower frequency ( fio/8), diminishing the frequency requirements for the utilized 


hardware components. 


The FEs are manufactured as a printed circuit board (PCB) that incorporates a 
Vivaldi antenna, a monolithic microwave integrated circuit (MMIC) and other 
electronic surface mount device (SMD) components and is housed in a metal 
case. The entire assembly measures only 11 mm x 35 mm x 50mm (W x H x 
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Figure 4.3: Photograph of receiver front end 


L). A photograph of the FE can be seen in Fig. 4.3. Brass metal was chosen 
for the housing material in order to mitigate the influence of surface oxidation 
which can decrease the electric conductivity of the case and thus weaken the 
ground connection between the FE PCB and the connected coaxial cable. The 
PCB was realized as a compound of four copper layers separated by dielectric 
substrates. For the top RF layer the halogen-free, ultra-low transmission loss 
material R-5515 from Panasonic with a relative permittivity of e, = 3.09 and a 
thickness of 127 um was used [Pan21]. With this, the width of the micro strip 
line could be set to 313 um to achieve a characteristic impedance of 50 Q. 


Antenna 


For the reception of the radar signal emitted by the RuT and the re-transmission 
of the modified signal, a co-planar Vivaldi antenna was designed [Bal15]. For 
the realization of the antenna, all layers of the PCB except the RF substrate 
were cropped, which can be seen in Fig. 4.3 by the slight glow in the antenna's 
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Figure 4.4: Simulated S-Parameters of the co-planar Vivaldi antenna 


surrounding area. The results of the S-Parameter simulation are shown in 
Fig. 4.4. The wide bandwidth of the antenna supports the applicability of the 
RTS system for the entire automotive frequency spectrum from 76 to 81 GHz. 
Its poor selectivity is tolerable considering its deployment in a controlled envi- 
ronment without other sources of electromagnetic radiation other than the RuT. 
An end-fire antenna type was chosen in order to reduce the electromagnetic 
backscatter of the geometric structure which may be interpreted by the RuT as 
a static target. 


RF Chip 


For the down- or up-conversion of the received or modified radar signal, a 
transceiver MMIC provided by the Leibniz Institute for High Performance Mi- 
croelectronics (IHP) was used [NHK19]. It incorporates both a receive as well 
as a transmit path including a low-noise amplifier (LNA) or power amplifier (PA) 
and an in-phase and quadrature (IQ) mixer and frequency multiplier each. As 
mentioned before, in order to mitigate mutual coupling and thus the generation 
of ghost targets, receive and transmit FEs are implemented as separate modules, 
which is why only the receive or transmit path of the MMIC is used. Moreover, 
the quadrature port of the IQ mixer was terminated and only a real-valued signal 
was transmitted in order to save hardware resources as will be further explained 
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Figure 4.5: Photograph of the radar transceiver MMIC 


in chapter 4.3. The chip was not originally designed for the deployment in a 
RTS system and could be adapted for this application. However, due to the 
efforts of radar chip manufacturers to integrate as many functionalities into a 
single MMIC as possible, the chip at hand was the best available choice at 
the time that provides direct access to the broadband intermediate frequency 
(IF) in- and output signals. The MMIC was placed in a cavity cutout of the 
PCB which allowed for relatively short bond wire connections. An enlarged 
photograph of the chip can be seen in Fig. 4.5. 


4.1.2 Back End 


For the digital back end of the RTS system, a ZCU111 evaluation board from 
Xilinx was employed, that features a Zynq UltraScale+ XCZU28DR radio 
frequency system on chip (RFSoC) with eight integrated high-speed ADCs and 
DACs [Xil18]. The Zynq SoC family from Xilinx integrates both a processing 
system (PS) and programmable logic (PL), in this case in the form of a quad-core 
Arm Cortex-A53, a dual-core Arm Cortex-R5, and an UltraScale+ FPGA. PS 
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Figure 4.6: Photograph of the Zynq UltraScale+ RFSoC ZCU111 evaluation board with mounted 
adapter board 


and PL are connected to their own dedicated Double Data Rate (DDR) Random- 
Access Memory (RAM) that is located on the evaluation board. Advanced 
eXtensible Interface (AXI) connections between the PS and PL enable the 
control of the FPGA unit through software applications that run on the ARM 
processors. A custom-made adapter board that is mounted on the ZCU111 and 
incorporates balanced-to-unbalanced (BALUN) converters enables the access 
of the ADCs and DACs through SubMiniature version A (SMA) connectors. A 
photograph of the setup is shown in Fig. 4.6. 
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Analog-to-Digital Converter (ADC) & Digital-to-Analog Converter (DAC) 


The eight integrated ADCs and DACs are configured to a common sampling 
frequency of fr = 4G Hz. The sampling clock is generated external to the 
RFSoC by an on-board RF phased-locked loop (PLL). The ADCs feature a 
resolution of 12 bit/sample and the DACs 14 bit/Sample, which are each scaled 
up to 16 bit/Sample for compatibility reasons. The maximum input respectively 
output power of the converters is Papc = Ppac = 1dBm. Since the sampling 
frequency exceeds the internal clock frequency ofthe FPGA to which the ADCs 
and DACs are connected to, the samples are parallelized at the digital output of 
the ADCs and serialized at the input of the DACs, using first-in-first-out (FIFO) 
buffers. The ADCs and DACs introduce an inevitable latency to the signal that 
passes through them, which together with other latencies of the system defines 
the minimum simulatable distance of a virtual radar target. Compared to others, 
the combined latency of ADC and DAC of tapc + Tpac = 162 ns is decisive 
and it alone increases the minimum target range by Rmin = 24.28 m. 


Field Programmable Gate Array (FPGA) 


The UltraScale+ FPGA of the XCZU28DR RFSoc features 850560 flip-flops 
(FFs), 425 280 lookup tables (LUTs), 1080 Block RAM (BRAM) blocks and 
4272 digital signal processor (DSP) slices. It receives the sample data stream 
containing the down-converted and digitalized radar signal from the ADC, ap- 
plies the target generating signal modifications and forwards it to the DAC for 
re-synthesis. The internal clock frequency of the FPGA is set to 500 MHz, 
which requires the data samples to be parallelized in groups of eight, to prop- 
agate through the PL simultaneously by means of an AXI stream interface. 
Consequently, each stage of the target generation unit on the FPGA must be 
implemented in parallel in order to avoid any backlog. 


Processing System (PS) 


The processing system hosts the ARM cores on which a baremetal C application 
is deployed. It controls the delay, frequency shift, attenuation and phase shift 
that the target generation unit on the FPGA applies to the radar signal. For 
this, it sets the corresponding registers in the PL through an AXI Lite interface. 
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In addition, the software application instances either a Transmission Control 
Protocol (TCP) or alternatively an User Datagram Protocol (UDP) webserver 
depending on the use case scenario. The server receives the list of targets that 
are to be simulated by the RTS system over an Ethernet connection from either 
a virtual environment simulation or user input. For this purpose, an open- 
source lightweight Internet Protocol (1wIP) stack is implemented which grants 
a resource-eflicient, low-latency transfer of the target generation parameters 
from the user or environment simulation to the FPGA. 


4.2 Dynamic Range 


The exact values of the transmit power and the receiver sensitivity of a com- 
mercial radar sensor are often times not published by the radar manufacturers. 
However, the absolute power levels of the receive and transmit chain of the RTS 
system, although not entirely negligible, are not as relevant in this application 
and can simply be adjusted with variable gain amplifiers and attenuators during 
the calibration process prior to the actual RTS measurements. It must only be 
ensured that the signal power at the receiver of the RuT does not drive its low- 
noise amplifier (LNA) into compression. More important in this context is the 
dynamic range of the system, as it is essential in order to simulate entire traffic 
scenarios that may consist of a large number of targets with various reflection 
properties. The system must be able to generate targets with a high RCS at 
small distances and targets with a low RCS at large distances at the same time. 
The signal power attenuation caused by a target at a distance of R with an RCS 
of o can be expressed as 


ra? 


CE (4.1) 


Atgt = 
The dynamic range that needs to be covered by the RTS system can be de- 
termined by calculating the difference of the extremes. For instance, the 
weakest radar reflection that is still detectable by a commercial automotive 


radar is caused either by a pedestrian (opa = —7 dBm”) at Rpa = 110m, a 
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moped (omp = 0 dBm’) at Rmp = 170m or a motorcycle (Ome = 5 dBm?) at 
Rme = 220 m [Con21], which all result in a signal attenuation of 


Argt,max|4p = -158 dB (4.2) 


The strongest receive signal that may impinge on the RuT originates from a 
high RCS target (truck, oj, = 20 dBm?) at the minimum simulatable range of 
Rmin = 25 m (cf. section 4.1.2) resulting in a signal attenuation of 


Atgt,minl yg = —105 dB (4.3) 
With this the required dynamic range can be determined to 
DRreq dB 7 Atgt,minlag = Atgt,max dB 7 53 dB (4.4) 


The integrated DAC features a resolution of 14 bit and therefore a dynamic 
range of 


DRpac|ag © 6.02 dB /bit - 14 bit = 84.28 dB (4.5) 


which sufficiently covers the required dynamic range of the system [PL15]. 


4.3 Spectral Model 


In the following, the spectral model of the developed RTS system is elaborated 
in order to facilitate the comprehensibility of the system design decisions made 
and to explain the individual processing steps within the RTS signal chain. 
The presented system was developed for current automotive radar sensors that 
occupy a bandwidth of B = 1 GHz at a carrier frequency of fer = 76.5 GHz 
[Con21, Rob21], but it can be easily adapted for future sensors with greater 
bandwidths or at different frequency bands, as only the sampling frequency of 
the ADCs and DACs or the LO frequency must be adjusted, respectively. The 
design of the RTS system is based on a double-heterodyne concept, meaning that 
the radar signal emitted by the RuT and received by the RTS receive antenna 
is down-converted in two steps; first analog to an intermediate frequency of 
fits = 1 GHz and then digitally to baseband. The up-conversion is done similar 
but in reverse. This two-stage process was chosen for two reasons. On the one 
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hand, the use of the intermediate frequency allows the signal to be digitalized 
and later synthesized in its real-valued representation, requiring only one ADC 
and DAC per RTS channel. On the other hand, the usage of the complex-valued 
baseband is required to generate virtual radar targets at a given velocity while 
suppressing ghost targets with the inverse velocity. 


A diagram showing the spectral progress of the radar signal through the individ- 
ual stages of the RTS system is shown in Fig. 4.7. Each plot shows the spectrum 
of the signal in GHz and is referenced by the block diagrams of both the front 
(Fig. 4.2) and back end (Fig. 4.8). The in-band amplitude slopes shown in the 
diagram indicate the directionality of the signal frequency band, which is of 
importance to the signal processing of the RuT. 


Fig. 4.7a shows the frequency band of the radar signal emitted by the RuT 
and received by the RTS FE at fer + B/2 = 77 GHz, and the LO signal that it is 
down-converted with at fio = 76 GHz, as a Dirac impulse. Both signals are real- 
valued, meaning their frequency components are mirrored in the negative half 
of the spectrum. The resulting signal (Fig. 4.7b) features frequency components 
at fits = 1 GHz and 153 GHz, however, the latter is not output by the MMIC. In 
order to mitigate aliasing, the signal is filtered using an analog low-pass filter 
(LPF) before it is digitalized by the back end ADC. Subsequently, the signal 
is digitally down-converted to baseband using a digital IQ mixer and a LO 
frequency of —1 GHz, and low-pass filtered using a digital LPF. The spectrum 
of the emerging complex-valued signal is shown in Fig. 4.7c, with indications 
for the spectral repetitions at + fs,rs/2 = +2 GHz. In the next step, the only 
spectrally relevant signal modification of the target generation process, that is 
the Doppler shift, is applied (Fig. 4.7d) and the signal is up-converted to its 
previous intermediate frequency of fits + fp = 1 GHz + fp. In order to save 
hardware resources, thereupon only the real part of the signal is continued, 
which brings back the mirrored frequency component in the negative half of 
the spectrum depicted in Fig. 4.7e. Next, the other target signal modifications 
are applied and the signal is re-analogized by the DAC and once more low-pass 
filtered to suppress the arising spectral repetitions (Fig. 4.7f). Finally, the signal 
is up-converted back to its original carrier frequency and re-transmitted towards 
the RuT by the transmitter FE. As can be seen in Fig. 4.7g, this results in an 
undesired spectral component at 75 GHz in addition to the intended band at 
for + B/2 = 77 GHz. However, this presents a tolerable side effect, since the 
unwanted frequency component, although also received by the RuT, is directed 
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Figure 4.7: Spectral progression of the radar signal through the RTS system with all frequencies 
given in GHz 
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in the opposite frequency direction over time ¢ and thus causes a de-chirped 
signal at the RuT that features a frequency slope double the steepness of the 
transmitted ramp 


B 
fo =2fus +27 +t = [2,4]GHz (4.6) 
S 


This signal component is filtered out by the LPF preceding the ADC of the 
RuT. Alternatively, the signal emerging from the target generation unit can also 
be re-analogized in its complex-valued representation using two instead of one 
DAC per RTS channel which avoids this behavior. 


The delay, attenuation and phase shift applied by the target generation unit do 
not have any effect on the spectral properties of the signal and can therefore be 
implemented as real-valued operations. The complex-valued baseband is only 
needed for the application of the Doppler shift, as otherwise an additional ghost 
target with the inverse velocity would arise. 


4.4 Target Generation 


The digital target generation unit is implemented using the configurable logic 
blocks in the FPGA. A block diagram displaying its functional components 
is shown in Fig. 4.8 with references to the spectral model in chapter 4.3. 
The signal, after being digitalized by the ADC, is down-converted to complex 
baseband using a digital down-converter (DDC). Subsequent to the application 
of the Doppler shift, which is realized with a direct digital synthesis (DDS) 
component and an IQ mixer, the signal is up-converted back to fis using a 
digital up-converter (DUC). After this, the delay, phase shift and attenuation are 
applied before the signal is re-synthesized by the DAC. As mentioned before, 
the sampling frequency of the ADCs and DACs exceeds the internal clock 
frequency of the FPGA with which the data samples are processed, therefore 
all of the signal modification functions had to be implemented to process the 
samples in parallel. If more than one target per channel shall be generated, 
the signal can be split after the DDC and re-combined just before the DAC 
minimizing the required hardware resources. 
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Figure 4.8: Block diagram of the back end target generation 


4.4.1 Digital Up-/Down-Conversion 


The digital up- and down-conversion of the signal prior and subsequent to the 
application of the Doppler shift, is done using an LO signal with a quarter of 
the sampling frequency fiosts = fs.rts/4 = 1 GHz. This allows the LO signal 
to be sampled at exactly four cyclic instances per period, which limits its value 
domain to 


Srts,lo = (1, j =], =j} (4.7) 
This in turn enables the output of the digital converters to be expressed as 
Y =X: Srts,lo = {XT +jxQ, —XQ +jxr, -xr — jxQ, XQ — jx} (4.8) 


where x describes the input signal. The real and imaginary parts of the samples 
are separated into two dedicated signal paths, thus the above function can be 
implemented by exchanging samples between the paths and flipping the sign 
bit of certain samples. A digital finite impulse response (FIR) filter succeeding 
the down-conversion is employed in order to suppress the undesired sideband. 
After the DUC, the imaginary signal path is terminated since for the following 
signal modifications the real-valued signal representation is sufficient and it 
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helps to save resources on the FPGA and potentially increasing the scalability 
of the target generation unit. 


4.4.2 Range Simulation 


The range simulation of the virtual target is realized by delaying the signal 
samples. This is done through sample buffering, meaning storing the incoming 
samples in on-chip BRAM blocks and reading them back a specified amount of 
clock cycles later. The blocks feature two read-write ports whose memory map 
addresses can be set individually allowing to write and read at the same time at 
arbitrary locations of the BRAM block. One of the ports is used only for write- 
and the other only for read-operations. Both ports cycle over the memory space 
with an address offset of 


Mrts = Trts * Isrts (4.9) 


between them which is determined by the desired signal delay tT, and the 
sampling frequency fs,rts of the ADC and DAC. A given range is simulated by 
calculating the required delay according to (3.1) and setting the correspond- 
ing address offset of the ports through the PS. The smallest realizable delay 
increment is defined by the sampling frequency 


1 
Ati = (4.10) 
rts Fas 
which in turn defines the range resolution of the RTS system to 
AR E aay (4.11) 
im = = á mm . 
sim Ds 


The maximum range that can be simulated is defined by the size of the memory 
space that is allocated for each RTS channel. 


4.4.3 Doppler Simulation 


The velocity of the virtual radar target is simulated through a complex-valued 
frequency shift. This is done by synthesizing an additional signal s;ıs,p With 
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the corresponding Doppler frequency and mixing it with the radar signal. The 
discrete complex-valued mixing product can be expressed as 


Y =X: Srts,D = XISI — XQSQ +] (x15Q + xası) (4.12) 


where x describes in the input signal. As mentioned before, complex-valued 
signals within the FPGA are realized as two designated signal paths (real and 
imaginary). DSP slices are employed for the fixed point multiplication and 
subsequent summation of the respective signal I- and Q-parts. The required 
Doppler signal is created using a DDS component with an adjustable phase 
increment that is accumulated in the phase register which is driven by areference 
clock. The register operates the address port of a LUT in which the values of 
a sin-function are stored. The output is the real-valued time-discrete sequence 
of a sin-wave whose frequency can be adjusted by setting the increment of the 
phase register. The imaginary part of the signal is obtained through a second 
output port of the DDS, whose address pointer is offset by a quarter of the 
wavelength. 


Given the velocity of the target, the required Doppler frequency fD,rts can by 
calculated according to (3.3) and the phase increment needed to synthesize the 
corresponding signal is 


Dart 


ref 


Aus = 27 (4.13) 
where fref describes the frequency of the reference clock. Since the smallest 
address increment is equal to one, the smallest phase increment 27 /NLur 
is defined only by the depth of the LUT Ni yr, as it must always span 27, 
whereas the maximum phase increment must not exceed z in order to meet the 
Nyquist-Shannon sampling theorem [Whi15]. Together with the reference clock 
frequency, the smallest phase increment determines the Doppler frequency 
rasterization 


= fref 


Nut 


A fo (4.14) 
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Given the depth of the LUT Niyr = 27° and the reference clock frequency 
fiet = 1 MHz, the velocity resolution of the RTS system can be derived as 


Co fref 


= 0.002 m/s 4.15 
2 ferNLUT / oy) 


AVsim = 


4.4.4 RCS Simulation 


The RCS of a virtual target is simulated by controlling the attenuation of 
the radar signal. This is achieved by multiplying the incoming samples with 
an adjustable factor that is stored in a register and set through the PS. For 
the multiplication, again a DSP slice is employed. The required amplitude 
attenuation is defined by the RTS channel gain that is needed to simulate a 
virtual target with a given RCS, which can be calculated according to (3.6). 
The attenuation can be derived as 


VO sj 
Arts = V Grts =. Acal (4.16) 


2 
R sim 


where Acai is the channel-specific amplitude offset that is determined during 
the calibration process. 


4.5 Resource Utilization 


The resource utilization is relevant, since it allows an estimation of the scala- 
bility of the design given the available resources of the utilized FPGA platform. 
Table 4.1 lists the number of logic blocks that are required for the implementa- 
tion of a single target generation unit and the RTS channel overhead structure. 
As stated before, not only does the back end feature multiple RTS channels, 
but also each channel may comprise of one ore more target generation units. 
Depending on the application it might be required to simulate multiple virtual 
targets whose corresponding radar echoes originate from similar angles, es- 
pecially when taking advantage of the concept of arbitrary angle simulation 
presented in chapter 6. Therefore, Table 4.1 distinguishes between resources 
needed for the generation of a single target and those common for the entire 
channel. The items AX/ and Combiner represent the configuration bus system 
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Table 4.1: FPGA Resource Utilization 


LUT FF BRAM DSP 


Channel 1313 4333 - 80 
AXI 190 141 - - 
LPF 738 3040 - 80 
Combiner 385 1152 - - 
Target 2843 3052 45 48 
AXI 152 139 - - 
Doppler 81 300 0.5 48 
Delay 458 651 4 - 


Atten. 2152 1962 - - 


and the merge unit that unifies the output signals subsequent to the individual 
target generation units, respectively. 


Table 4.2 enumerates the resources available on the Zynq UltraScale+ 
XCZU28DR along with those necessary for the basic system functionality 
(Fixed), such as PS-PL interconnects and the system reset control. In addition, 
the number of logic blocks required to realize eight RTS channels and a total 
of 32 target generation units are listed. The number of channels is chosen with 
respect to the number of ADCs and DACs integrated on the chip. The number 
of target generation units is selected exemplarily to demonstrate a potential use 
case, where each channel may simulate up to four virtual targets simultaneously. 


As can be seen, with the given parameterization, the proportionate usage lies in 
the low tens for all logic block types, except for the DSP, which is employed for 
the multiplication operations of the FIR low-pass filter (LPF) and the complex 
mixing of the Doppler unit. Although the DSP demand can be reduced, this 
comes at the cost of a drastic increase of the LUT and FF usage [dDP09], hence 
not necessarily more targets can be simulated. Moreover, it is not advisable to 
scale the design beyond a usage rate of 50 % since all radar-signal-related data 
paths are clocked at a high frequency close to the FPGA’s maximum, making 
their routing time-critical and complex [X1122]. 
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Table 4.2: FPGA Design Scalability 


LUT FF BRAM DSP 


Available 425280 850560 1080 4272 


Fixed 6222 4951 - - 
8x channel 10504 34664 - 640 
32x target 58976 97664 144 1536 
Total 75702 137279 144 2176 
Usage 17.8% 16.14% 13.33% 50.94 Yo 


4.6 Signal Model 


In the following, a generic signal model that analyzes the propagation of the 
radar signal from the RuT through the RTS system and back is presented. It 
further depicts the successive radar signal processing within the RuT and dis- 
closes the impact of the various RTS parameters on the target detection. The 
model is based on the analytical derivations of a chirp sequence radar in section 
2.2.1, but is adapted in order to examine the influence of certain RTS charac- 
teristics, which in turn allows to deduce the target generation configurations 
required, given the properties of the real-world radar target to be simulated. As 
before, the model is represented mainly by its phase, with its amplitude mostly 
neglected. It will later be re-used and adapted for the application-specific use 
cases in chapter 5 and 6. 


The signal emitted by the RuT and its phase can be expressed as 


+ O(n) , t € [0, Ts] (4.17) 


B 
tne) =2 epa «e 
Pix (t, Ne) T A IT, 


Six(t, Nc) = exp {jpx(t,n.)} = exp {ian fe “t+ x Pl + on) (4.18) 
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It travels from the RuT through free space to the RTS FE, where it is received 

with a time delay of 

R 

= (4.19) 
0 


Tx = 


where R. describes the physical distance between the RuT and the FEs, which 
is assumend to be the same for the transmit and receive paths. Subsequently, 
it is down-converted using an LO signal with a frequency of fio, whose phase 
can be described by 


Plo(t, Nc) = 27 fio (t + ncTe) (4.20) 


to an intermediate frequency of 


Sits = Ser = fio (4.21) 


Thus, the signal phase at the output of the RTS FE can be expressed as 
Prts,rx(t, Ne) = P(t — Tix, ne) — Ylo(t, Nc) (4.22) 


B 
= Qn — forTix + Tits “f+ ay = Tx)? = Fo e iele + D(n¿) 
S 


Next, the signal is digitalized using the integrated ADC of the RFSoC. There- 
fore, the time domain is discretized using the sampling frequency £ rts 
Mrts 
t> —, Ms € No (4.23) 
s,rts 
where nys denotes the sample index, resulting in a signal phase at the output of 
the ADC of 


2 
Arts B Arts 
Prts ade (Arts N.) = 2n |—fortix tJe —— + = | = ra) 
? Fer 2T; fons 


— fio NeT¿| + ©(n.) (4.24) 
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Thereupon, the Doppler shift is applied to the signal. For this, the DDS unit 
synthesizes a signal with a frequency of 


2Vsim 


A 


JD.rts = (4.25) 


that corresponds to the radial velocity of the simulated target. The synthesized 
signal with a phase of 


n 
Prts,dds (Arts» ne) = 27 fp, rs (>= + ner) (4.26) 


s,rts 


is mixed with the radar signal and the resulting signal phase at the output of the 
Doppler unit can be expressed as 


Prts,D (Mets, Nc) = Pris,ade (Arts ne) = Prts dds (Arts» Nc) (4.27) 


=2n =r + (fits — fo, rts) 


Nts 


Ss, rts 


B : 
| = + O(n) 


+ = 
2T; Fons 


2 
= ra) = (fio + JD.rts)ncTe 
After that, the signal is fed into the delay unit, where it is buffered for 7,5 
samples, which corresponds to a delay of 


2Rsim _ Thts 
<= ale S 
co frits 


Mrts € No (4.28) 


Trts = 


The signal phase subsequent to the delay unit can be described by 


Pris,r (Arts, Me) = Prts,D (Arts — Nets» Ne) (4.29) 
Mrts 
= 27 | fertix — (fits — fd,rts) Ts + (Lets — JD.rts) 
Fsrıs 
B jn 2 
ts 
“IR (7 i tus) - (fio + forts) McTe + O(n) 


where the number of buffered samples 7,1, is substituted by the corresponding 
delay T,1s. As the final step of the target generation unit, the RCS is simulated. 
For this, the signal is attenuated by multiplying the signal samples with a factor 
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that can be calculated according to (4.16) and that corresponds to the range and 
RCS of the target 


Srts,A (Arts, Nc) = Arts * Sris,7 (Arts, ne) (4.30) 


Subsequently, the signal is converted back to the analog domain, re-substituting 
the discrete sample index with a continuous time variable 


Mrs > Ís.1ts -t (4.31) 


The signal phase at the output of the DAC can be expressed as 
Prts,dac (t, Nc) =21 | = FerTix = (Arts = Jp, rts) Tets + (fits = Jos) -t (4.32) 


B 
+ I. (t — Tx — Tr) = (fio $ Jors) -ncTe| + D(n¿) 
s 


Next, the signal is fed into the transmitter FE where it is up-converted back to its 
original carrier frequency using the same LO signal as before and re-transmitted 
towards the RuT. The signal propagates through free space and is received by 
the RuT with a delay of 


Ti == (4.33) 


Combining the delay caused by the free space propagation to and from the RTS 
into a physical delay of 


2R 
To = Tix + Tx = — (4.34) 
co 


and this in turn together with the artificial delay of the RTS into a common 
delay of 


T = Te + Tits (4.35) 


64 


4.6 Signal Model 


allows to simplify the resulting signal phase at the receiver of the RuT to 
Prx(t, nc) = Prts,dac (t = Trx, Me) + Plo(t — Tix, Nc) (4.36) 


=2n| - FerTe = (fits > JD, 1ts) Tets + JD, rtsTrx + (Jor = fp,rts)t 


B 
+ 37," - T)? — fDirts © NcTe| + ©(n.) 


The beat signal is formed by down-converting the receive signal with the trans- 
mitted one which translates to a subtraction of their respective phases 


volt, Ne) = Pix(t, Ne) — Prx(t, Ne) (4.37) 
B 
=2n (Fr = fou] EF JD.rts "Nele 
s 


B 
+ FerTe + (Fis = fo,rts) Tits — Jd, rtsTrx = a 
s 


4.6.1 Radar Signal Processing 


In the following, the radar processing of the signal returning from the RTS is 
elaborated. As in chapter 2.3, the expression for the phase of the beat signal is 
simplified using the relations in (2.17), (2.18) and (2.19) to 


B 
p(t, Ne) = 27 T oe Jd, rts "Nele + FerTe + SitsTrts (4.38) 
s 


Range-Doppler Estimation 


Passing the signal described above through the same processing stages as in 
section 2.3.1 and 2.3.2 allows to express the output of the two-dimensional 
range-Doppler FFT as 


sp(k,m) = NsN. - sinc (Bt — k) - sine (fp,rts Im — m) (4.39) 
. Br-k  forsIim-m 
- exp fiar — + a an + forT, + fatas} 


65 


4 Test System Setup 


From this, the respective range- and Doppler cell that holds the target can be 
identified as 


Kdet = Br (4.40) 
Mdet = FD.rtsTm (4.41) 


which in turn enables the estimation of the range and velocity of the virtual 
target to 


co 
Raet = Kaet h 2B = Re + Rsim (4.42) 
A D, rts 
va =m 5 = POA Z vm (4.43) 


As can be concluded from (4.42), Rsim must be adjusted to account for the 
physical distance between the RTS and RuT in order to simulate the correct 
target distance. In addition, Rsim must be calibrated for each RTS channel 
individually, since their respective signal path lengths within the FPGA logic 
may vary from another which in turn introduces differing latencies. 


Angle of Arrival Estimation 


Taking into account the lateral position of the FE from where the radar signal 
is re-transmitted towards the RuT from, the free space signal delay can be 
redefined as 


Re  dısin(0) 
— + ———— NA 


co co 


Trx = 


(4.44) 


where 6 € [-90°, 90°] describes the azimuth angle of the transmit FE from 
the perspective of the RuT. The signal phase subsequent to the range-Doppler 
detection can be expressed as 


da sin(@) 
a 


B B 
gp (na) =2n (fs + 3) To + ES + 4 Tris + NA (4.45) 
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4.6 Signal Model 


Applying Fourier beamforming and simplifying in the same manner as in section 
2.4.2 results in a signal of 


sala) = Na - sinc 4 [sin(9) — sin(a)] ma) (4.46) 
"ESP {i2n B $ 3) To + Im + 3) Trts + faso) ae || 


from which the angle of the RTS FE and thus the AoA of the virtual target can 
be estimated 


sin (ae) = sin (0) (4.47) 
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5  High-Precision Range Simulation 


The necessity of simulating virtual radar targets at a high-precision range is not 
intuitively apparent. Common RTS systems are already capable of adjusting 
the range of virtual targets in increments in the single-digit centimeter region 
[SN20,GGS*17], which seems sufficient to deceive current radar senors with 
range resolving abilities of AR = 15cm [Con21,Rob21]. Moreover, given the 
context of automotive radar target simulation, it misleadingly appears irrelevant 
if a preceding vehicle is simulated with a rage offset of a few centimeters. 


However, not only does the range accuracy of a modern sensor exceed its 
resolution and therefore must be covered by the validating system, but also 
will future radar sensors utilize the entire automotive radar frequency spectrum 
between 77 and 81 GHz [Eur08], which implies an increase of the radar signal 
bandwidth and thus a finer range resolution of AR = 38 mm. Moreover, real- 
world objects, e.g. a car, usually cause a multitude of reflections, forming so 
called point clouds that consist of many individual radar targets close to each 
other in the lateral and longitudinal domain. Therefore, RTS systems must be 
capable of generating multiple virtual radar targets with small radial offsets 
which in turn requires a fine range precision. Furthermore, fast-moving objects 
feature characteristics that make them too complex to be substituted by ideal 
virtual point targets as generated by the classical RTS approach. These objects 
cause a migration effect in the signal processing of a radar sensor that degrades 
the target detection and must be replicated for credible radar target simulation. 
This necessitates the capability to adjust the simulated target range seamlessly 
even within a measurement cycle. 


The aforementioned challenges can be overcome with an improved range sim- 
ulation precision. The precision at which virtual radar targets can be created 
by an RTS system in the range domain is driven by the rasterization at which 
the delay that the system applies can be adjusted. It is obvious that the smallest 
possible delay step limits the range simulation precision and must therefore be 
minimized. It also determines the update period at which the range of a virtual 


69 


5 High-Precision Range Simulation 


target has to be re-set given its radial velocity. Updating the simulated target 
range within a measurement cycle at a high refresh rate allows to replicate the 
range migration phenomenon and enables the generation of credible virtual 
radar targets. The adjustment of the Doppler shift, that represents the velocity 
of the virtual target, does not pose an overly difficult challenge, since target 
velocities do not change as quickly as ranges and the components utilized for 
this task already possess fine rasterizations. Moreover, the synthesis of range 
migration inevitably leads to the same effect in the velocity domain, as will 
later be proven [DAN*21]. 


In the following, first the necessity of high-precision range simulation is sub- 
stantiated which enables the synthesis of the migration phenomenon that can be 
observed in the presence of high-speed radar targets which are common in the 
automotive context. Next, the generic signal model in chapter 2.2.1 is extended 
in order to facilitate the comprehension of the migration phenomenon, analyze 
its origin and estimate its magnitude. Thereupon, a reference measurement 
conducted with a magnetic sledge is presented that verifies the preceding an- 
alytical elaborations. After that, the range increment and update rate required 
to synthesize the migration effect in radar target simulation is examined and 
realized using fraction delay filter. Subsequently, a verification measurement is 
shown that validates the successful implementation of the migration synthesis 
in the RTS system. Finally, some practical considerations are given. 


5.1 Migration Phenomenon 


Conventional range-Doppler processing with a coherent processing interval as 
described in chapter 2.3 implicitly assumes that the range change of a dynamic 
target within one measurement period due to its radial speed is less than the 
range resolution. The same is adopted for the velocity estimation, that is subject 
to a spread caused by the dependence of the Doppler shift on the instantaneous 
frequency. The migration phenomenon describes the spread of a single radar 
target over several range and Doppler bins. It is primarily known from synthetic- 
aperture radar (SAR) applications [PDF99], but is also gaining importance in the 
automotive sector, due to the continuous increase in automotive radar resolution 
capabilities, which pronounces the effect [XBP19]. 
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5.1 Migration Phenomenon 


Dari 


ÒR = vIm 


Figure 5.1: Diagram of the range migration phenomenon principle in the context of an automotive 
application 


Fig. 5.1 shows exemplarily the principle of range migration, where the radar 
target travels the distance ôR in the measurement time Tm. The drawing is by no 
means to scale and is for illustrative purposes only. As can be seen, the distance 
of the vehicle to the RuT increases during the measurement time, and so do 
the signal delay and the detected target range that is derived from each chirp. 
If this increase in the range R exceeds the range resolution capability of the 
radar sensor AR, the target’s peak in the range-FFT spreads over several bins 
and range migration can be observed. In order for this to happen, the product of 
the target’s velocity and the measurement time must be greater than the RuT’s 
range resolution 


co 
ôR = vIm > AR = — 5.1 
Vim 2B (5.1) 


To prove the occurrence of the migration phenomenon in automotive radar ap- 
plications and therefore the necessity to synthesize it for radar target simulation, 
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the minimum velocity, that is needed for the target range to migrate over a single 
bin, can be calculated. This can be done exemplarily for acommon automo- 
tive radar sensor, whose measurement time can be derived from its velocity 
resolution Av 


co 


= 
NTE 


(5.2) 


Considering a velocity resolution of Av =0.1 m/s, a carrier frequency of 
fer = 77 GHz, a radar bandwidth of B = 1 GHz [Rob21] and substituting Tm 
in (5.1) reveals the needed velocity of 


co = Le, = 7.7 m/s 7X 28 km/h (5.3) 


> = 
“* ORT, B 


It is understood that velocities in this order of magnitude occur on a regular bases 
in road traffic and thus it can be concluded that the migration phenomenon plays 
a role for the validation of automotive radar senors and the ADAS functions 
that depend on them. 


5.2 Extended Chirp Sequence Signal Model 


In order to facilitate the comprehension of the migration phenomenon and to 
unveil its analytical origin, a detailed signal model is developed. As before, the 
signal is represented mainly by its phase, as this neatly illustrates the various 
effects of the radar target. The model is based on the derivations in chapter 2.2.1, 
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but leaves out some of the simplifications made there. Considering (2.15), the 
beat signal phase can be expressed as 


Range-Doppler coupling 


— 
Yp(t, Nc) = 2 ($e (1 _ 22) + 2 -t+ (2 (Fe _ 27) "Nele 


U 
U 
Range estimation Velocity estimation 


B2 B 2y? 
t ri 1-2))r- T s - (neTe)” 
T; co co Ts c 


+(#2 (1- 22) t+ NcTe + To (da e) (5.4) 
A ee 


Range-Doppler migration Angle of arrival estimation 


The expression can be simplified with regard to the range-Doppler migration 
and the range and velocity estimation, leaving out the range-Doppler coupling 
and the constant phase term. The terms that depend on 1? and ne? are also 
negligible, since they present a much smaller impact on the signal phase than 
the aforementioned ones. In addition, the relations in (2.17), (2.18) and (2.19) 
are considered to streamline the expression even further. After the discretization 
by the RuT’s ADC, the phase can be formulated as 


B 2 B2 
yo -ns + aid "Nelc+ er ns: ne (5.5) 
Ss 


Pb(Ns, Nc) = 27 7 Na 


Applying the range-FFT and then simplifying in the same manner as in section 
2.3.1 leads to 


2 
A ~ 4) (5.6) 
co 


2 
- exp 3327 T "Nele + 


da -) 
2 


As can be observed, in this equation the argument of the sinc-function also 
holds a delay offset, that is caused by the range change of the moving target 
(Rmig = vn.Te) and that increases from chirp to chirp. The detection range bin 
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kaet is subject to the influence of the chirp index n., and therefore the detected 
range of the last chirp differs from that of the first chirp by ôR 


2v 
kar = B [m u ne) (5.7) 
co 
R =v; (Ne — DT. = vTm (5.8) 
The number of bins that the range migrates over is given by 


ôR 2BTm 
= v 


Naes AR a 


(5.9) 


This number has to be considered with +1, as the number of bins visited also 
depends on the position of the initial range peak within the bin. 


The Doppler-FFT is applied directly to (5.5) and not subsequent to the range 
processing, since the FFTs can be exchanged due to the distributive law, and it 
is more demonstrative to show the migration effects in this representation. 


2v B 
,m) = Ne- si — Tm |1 + ==: — 5.10 
sp(ns, m) a sine 7 a + Na n) m) ( ) 
To (1+ NEE ns) - m 
exp 3327 | — To ` ns + > 


This time, it can be witnessed that the maximum of the sinc-function depends 
on the sample index ns and thus the detected velocity experiences a spread over 
multiple bins 


2v B 
= —T,,|1 . Al 
Meet qm | T Nofa n) 65.11) 
B 
=v{l+——.- 5.12 
Vdet = V | N, fa n) ( ) 
with a span of 
B B 
ôv = N; — 1) v x —v (5.13) 
Ns fer ( ) Je 
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5.2 Extended Chirp Sequence Signal Model 


Table 5.1: Simulation radar parameters 


Name Symbol Value 
Bandwidth B 1GHz 
Sampling rate fs  40MS/s 
Chirp repetition period Tr 30 us 
Number of chirps Ne 1024 


The number of velocity bins visited equals that of the range migration in (5.9) 


Nmig = 7, = v= (5.14) 


To illustrate the phenomenon and its effects, a numeric simulation was per- 
formed with the parameterization of a modern radar sensor [NXP20] shown in 
Table 5.1. For this, a target was simulated at an initial distance of R = 30m 
and with a radial velocity of v = 80 km/h = 22.2 m/s. The results of the simu- 


Norm. ampl. in dB 
Norm. ampl. in dB 


Range in m Velocity in m/s 


(a) Range profile (b) Doppler profile 


Figure 5.2: Simulation results of the respective one-dimensional (a) range FFT for all N. chirps 
overlaid and (b) Doppler FFT for all N, samples overlaid 
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lated signal after the respective one-dimensional FFTs can be seen in Fig. 5.2. 
In the range respectively Doppler profile, each trace corresponds to a chirp 
respectively a sample. The number of bins that the target migrates over can 
be calculated to Nmig = 4.55, which is consistent with the range and Doppler 
profile plots. 


5.3 Reference Measurement 


In order to validate the effect of the migration phenomenon, a reference mea- 
surement with amoving radar target was conducted and published in [DAN*21]. 


MN Rul AWR1443 
JE 


LEN 
De 


Figure 5.3: Photograph of migration reference measurement setup © 2021 IEEE [DAN*21] 
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5.3 Reference Measurement 


For this, a corner reflector was mounted onto a sledge that moves along a linear 
magnetic rail and that was accelerated to v = 1.25m/s. As the RuT, a Texas 
Instruments AWR1443BOOST radar sensor evaluation module was employed, 
which transmitted a chirp sequence signal with a bandwidth of B = 1 GHz ata 
carrier frequency of fer = 77 GHz. Since the maximum velocity of the sledge 
and the number of chirps were limited, the chirp repetition period was extended 
beyond what is typical for an automotive sensor to Te = 1.25 ms, in order to 
increase the measurement time Tm, and therefore to be able to showcase the 
migration phenomenon. Since the velocity of the sledge exceeded the maxi- 
mum velocity resulting from the parameterization chosen (Vmax = TT) the 
unambiguous interval of the velocity estimation was shifted to [Vmax, 3Vmax]. 
The number of chirps was set to N. = 240. A photograph of the measurement 
setup is displayed in Fig. 5.3. 


The results of the measurement are depicted in Fig. 5.4. Again, each trace 
corresponds to a chirp respectively a sample in the range respectively Doppler 
profile. It can be clearly seen that in both the range and Doppler profile the 
target migrates over several bins. The additional peaks in the range plot can 
be attributed to the metallic frame of the sledge and the cable routing that 
cause unwanted reflections. According to (5.9), the expected number of bins 


= = 
E # 
E = 
E 3 
E E 
9 9 
Z Z 
12 18 125 18 13 
Range in m Velocity in m/s 
(a) Range profile (b) Doppler profile 


Figure 5.4: Measurement results of the respective one-dimensional (a) range FFT for all N. chirps 
overlaid and (b) Doppler FFT for all N, samples overlaid © 2021 IEEE [DAN*21] 
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over which the target travels can be calculated to Nmig = 2.5, which is in 
agreement with the measurement results, considering Nmig has to be adjusted 
by +1 depending on the position of the initial range and velocity peak within 
the respective occupied bin. 


5.4 High-Precision Radar Target Range 
Simulation 


From the findings above, it can be concluded that RTS systems are required 
to create virtual radar targets with a high range accuracy. In the following, an 
extended signal model that is based on that of a generic RTS in chapter 4.6 is 
developed. It demonstrates how synthetic radar targets can be generated with a 
fine discretization in the range domain and allows an investigation of the update 
rate needed at which the target range has to be adjusted. As will later be shown, 
the delay increment that defines the smallest possible range offset undercuts the 
capabilities of sample buffering, for which reason fractional delay filters are 
introduced. Finally, the results of a measurement campaign are presented that 
verify the practical functionality of the proposed scheme. 


5.4.1 Signal Model 


Referring to (4.29), the artificial delay applied by the RTS back end is re-defined 
as 


Trts = Trts,0 + Tu (5.15) 
to consist of a constant part 
Mrts,O 
Trts,0 = a Nrts,0 € No (5.16) 
Fsrts 


with the index of the sample delay n.ı,o and a dynamic one that represents the 
target’s range change during the measurement 


2Ru _ 2Vsim 
co co 


“Nay (5.17) 


Tj = 
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where 7, describes the range update period and ny € [0, Nu — 1] the update 
index. This leads to a sample delay of 


Mrts = Trts Ls rts = (Trıs,o + Tu) Fsrts > Mrts € Q* 


2Vsi 
= [mu + E è m| Ssrts (5.18) 


Therefore, the resulting phase of the radar signal after the artificial delay has 
been applied can be expressed as 


Orts, 7 (Arts, Mc, Ny) = Prts,D (Mets — Nrtss Nc) (5.19) 


2Vsi 
— ferTix = (fits =, Jome) [mu + ~ i nl 
0 


2 

B Mrts 2Vsim 

Er — Tx — | Trts,0 + “nalu 
Ss, rts co 


+ (fits — Po, oe 


= 21 


— (fio + Sports) McLe + D(nc) 


am 


The subsequent conversion back to the analog domain and up to its original 
carrier frequency and the mixing with the transmit signal of the radar leads to 
a beat signal phase of 


B 
Pb(Ns, Nc, Mu) =2n (2»- fo, az F F fo, rts Nele 


i 
B 2Vsi 
au Im 2 JD.rts = ai) ==. Aylu 
Is co 
B 2vsim : Ms E B Am -(n aTa)? 
Ts co Ís Ty 


B 
+forTe + (Sets = Ports) Trts 7 JD.rtsTrx = a (5.20) 


S 
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Taking into account the relations in (2.17), (2.18) and (2.19) and To = Te + Trts,0 
the expression can be simplified to 


B Ng 2Vsim 
Po (Ns, Nc, Nu) =m | m t JD.rts "Nele + Jis — “Mula 
T fs co 
B2Vsim Ns 
+— — -— -n 7, (5.21) 
TG m 


After the application of the range-FFT and the same simplifications used before, 
the signal can be represented by 


2 el 
sed, = N sine 8 + “nT e) (5.22) 
co 
o [om p Ble 2B mK 
. EX — — 
Pp3J47 Netc 2 N, 


From this, the detected range bin, and therefore the estimated range can be 
expressed as 


2 A 
kat = B (v pom, mi (5.23) 
Co 
co 
Raet = Kaet * 2B = Ro + Vsim * Mula (5.24) 


The detected range comprises of a spread Ry = Vsim * 471 Which is a function 
of the update period and index. In the following, investigations are made on 
how often or seldom the simulated target range has to be adjusted, and therefore 
how T, has to be chosen, in order to synthesize the migration phenomenon in 
radar target simulation. 


5.4.2 Range Update Rate Investigation 


It is evident that the target range, respectively the delay, should be updated as 
seldom as possible in order to safe resources in the target generation unit. More 
frequent updates require greater memory allocation, as will be shown later. 
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Still, the same range change ôR as that of a real target must be covered during 
a measurement 


ÔR = Nulu = NeT (5.25) 


From (5.24) it follows thatthe simulated target range must be updated frequently 
enough for every bin to be visited that the target migrates over. For this, the 
range increment of a single update must be less than or equal to the range bin 
width 


Vsim * Ty < AR (5.26) 


Considering the maximum unambiguously detectable velocity of a target Vmax, 
the following can be concluded 


AR 
= Fer, (5.27) 


< 
di Vmax B 

Setting the update period to five or ten times the chirp repetition period 
(Ty = 5T. VT = 10T.) as an example and thereby fulfilling this criterion yields 
the desired migration phenomenon of the virtual target in the range and Doppler 
domain, as can be seen in Fig. 5.5a and Fig. 5.5b. This time the profiles only 
show the maximum values of each FFT bin for all samples or chirps respectively 
with a marker to indicate the actual bin value to allow an easier comparison. 
However, choosing the update period like that also leads to spurious peaks in 
the Doppler profile, as can be observed in Fig. 5.5c. The number of unwanted 
peaks correlates to the factor by which the update period exceeds the chirp 
repetition time. The ideal signal according to the signal model from 5.2 is 
plotted as a reference. 


These erroneous results can be illustratively explained by plotting the respective 
beat signal value and phase as a function of the chirp samples n. and for a 
single time sample n, as shown in Fig. 5.6. As can be seen, the beat signal 
and its phase stay constant between two range updates for the update periods 
of T, € (5T¿, 10T.), which creates harmonic peaks in the frequency analysis 
comparable to the spectrum of a digital signal. Again, the ideal signal according 
to the signal model from section 5.2 is plotted as a reference. 
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Norm. ampl. in dB 
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Figure 5.5: Range and Doppler profile simulation results for different range update periods Ta 


Reducing the range update period to equal the chirp repetition time prevents 
any steps in the beat signal phase and thus undesired, spurious peaks in the 
Doppler profile. Also, from the analytical perspective this makes sense, since 
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Figure 5.6: Simulation results of the beat signal (a) phase and (b) value for different range update 
periods Ty 


determining Ta = T, allows to streamline the expression for the detected range 
and velocity to 


Raet = Ro + Vsim * NeTe (5.28) 


B Sits B 
Vdet = Vsim | 1 + == 


fer Ng Ser l 
which largely correspond to those of the extended signal model in section 5.2. 
The detected velocity deviates slightly, however, this can be accounted for when 
setting Vsim- 


Ns (5.29) 


For these reasons, it can be concluded that the simulated range and thus the 
applied artificial delay must be adjusted at least once per chirp. However, 
the update period is not limited downwards and can also be set smaller than 
the chirprate (Ta < Te). The chirp repetition time can be derived from the 
maximum unambiguous velocity Vmax, as this parameter is usually given in the 
data sheet of a commercial radar sensor. For practical reasons, it is advisable 
to set the update period to half of the chirp repetition time (Ta = 7./2) in order 
to guarantee the suppression of any spurious Doppler peaks, while still limiting 
the resource demand. 
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5.4.3 Delay Increment 
The delay increment that is needed to adjust the target range after every chirp 
(Ta = T,) can be expressed as 


QAR, DIV 
co u co 


Ans (5.30) 


Considering the necessity to synthesize range-Doppler migration for all de- 
tectable velocities, the smallest delay increment required can be calculated with 


Vsim = Vmax tO 

21 -Vmax _ 1 
co 2 fer 

The smallest possible delay increment in digital RTS systems, which implement 


sample buffering for the range simulation, results from the sample frequency 
of the ADCs and is too coarse for delay steps of such fineness: 


(5.31) 


ATu,min = 


1 
ATs = — > ATu,min (5.32) 
Ss,tts 


Therefore, the target generation is extended with fractional delay filters that 
allow the application of artificial delays, even of fractions of the sample interval. 


Fractional Delay Filter 


The underlying theorem for applying a delay of a fraction of the sampling pe- 
riod to a discrete-time signal is based on the principle that a bandlimited signal 
can be reconstructed exactly if it has been sampled correctly. In this context, 
this implies that the sample frequency must be greater than the bandwidth of 
the signal (fists > B) [Sha84]. The idea is to convert a time-discrete signal 
Sin(Myts) to its time-continuous representation (x(f)), apply an arbitrary delay 
Ttg and convert it back to the time-discrete domain (Sout(Mrts)). The recon- 
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struction operation of the time-continuous signal can be expressed with the 
Whittaker-Shannon interpolation formula 


= t-T, n 

ý — F£s,rts * /rts 

x(t) = B Sin (Arts) * SINC (pee) 
Ts rts 


Mis 2700 


(9.33) 


where T, rts denotes the RTS’s sample interval. The continuous signal y(r) is 
subsequently time-shifted with a fractional delay Tja 


oo 

; t — Ta — Ts rts * Arts 

y(t) =x(t— Ta) = Y sins) sinc (os 
Ts rts 


Mrs 200 


) (5.34) 


Since the time shift is applied in the continuous-time domain, the delay r;a can 
be adjusted steplessly. Thereafter, the signal is discretized again at equidistant 
time instances t = Ts rts + Mts 


oo 


Sout (Mrs) = > Sin(Mrts) . sine 


Arts 00 


To rts * Mrts — Tfd Z T; rts * Arts 
Ts rts 


) (5:39) 


The expression can be simplified by setting 7; rs = 1 s without loss of generality, 
since the actual sampling rate does not need to be fixed for a discrete-time system 
description [DCD11], resulting in 


oo 


Sout(Mrts) = > Sin (Arts) - sinc (Mus — Tfd,d — Nrs) (5.36) 


Mrts=—0O 


where Tfaa = Tra/s describes the unit-free fractional delay. Comparing the 
equation to the definition of the discrete convolution, it becomes apparent that 
the operation can be implemented as a filter. However, the ideal fractional delay 
operation can not be realized since its impulse response is of infinite length. 
Therefore, an approximation needs to be made by using a finite impulse response 
(FIR) filter limiting the number of taps Nap by windowing the coefficients. 
Additionally, the coefficients are shifted by (Niap — 1)/2 to make the filter 
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Figure 5.7: Simulated exemplary calculation of the fractional delay filter coefficients © 2021 IEEE 
[DANZ21] 


causal, which in turn introduces additional latency to the signal. The coefficients 
can be calculated using 


> (5:37) 


h(nap) = sine Mtap — Tfd,d Z 


where Map € [1, Niap] is the filter tap index. Fig. 5.7 shows as an example the 
filter coefficients for no delay and Tja = 0.3 Ts ts with Niap = 9 taps. As can be 
deduced from Fig. 5.7a, the filter introduces an additional, inherent delay of 
(Map — 1)/2 = AT, ts, that deteriorates the minimum simulatable target range 
of the system. Prominently, the delay grows with the amount of filter taps 
employed, which must therefore be kept minimal. 


However, reducing the number of taps causes ripples in the filter transfer func- 
tion which in turn leads to a non-uniform frequency response that can inject 
ghost targets in the radar signal processing. Determining the number of taps of 
the fractional delay filter is therefore a compromise between these two consider- 
ations. Fig. 5.8 shows the frequency response of the filter for different number 
of taps Nap and for a set delay of Tia = 0.57; rts, which causes the strongest 
fluctuations in the transfer function. The filter is designed as an all-pass fil- 
ter, i.e. the bandwidth of the radar signal is smaller than the filter’s cut-off 
frequency and thus only the transfer function and its ripples that lay within 
the occupied signal band are of relevance. The influence of the fluctuations 
can be reduced by windowing the coefficients, as is exemplarily shown in the 
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Figure 5.8: Simulated frequency response of the fractional delay filter for different number of taps 
Map without (solid) and with (dashed) Blackman window 


graph for a Blackman window. Furthermore, the cut-off frequency itself is also 


dependent on the number of taps used, which has to be taken into account when 
designing the filter. 


5.4.4 Measurement Results 


In order to verify the functioning of the proposed approach to synthesize the 
range Doppler migration phenomenon in the context of radar target simulation, 
a measurement campaign was conducted and published in [DANZ21]. For this, 
the RTS system that is described in chapter 4 was set up in front of a Texas 
Instruments AWR2243BOOST radar sensor, that was utilized as the RuT, at a 
physical distance of Re = Im. A photograph of the setup of the measurement 
can be seen in Fig. 5.9. The receive and transmit front ends of the RTS 
were placed vertically aligned and behind a metal sheet that was used to block 
off undesired radar reflections of any structures in the background. The carrier 


frequency of the RTS was set to fits = 500 MHz. The sensor was parameterized 
according to Table 5.2. 


First, a reference measurement was performed for which an artificial radar 
target was generated with a range of Rsim = 40m and a radial velocity of 
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Figure 5.9: Photograph of the migration measurement setup 


Vsim = 25 m/s. No fractional delay unit was employed at this time and the basic 
back end design presented in chapter 4.1.2 was used. The range and Doppler 
profiles are shown in Fig. 5.10, where it can be seen, that the virtual target 
occupies only a single bin in each plot. The range of the target is subject to a 
small offset which originates from an imperfect range calibration. 


For the subsequent verification measurement, the same virtual target parameters 
were used, but the fractional delay unit was put in place. For this, a finite impulse 


Table 5.2: Radar parameters 


Name Symbol Value 


Carrier frequency fer 77 GHz 
Bandwidth B 1GHz 
Sampling rate fs 25MS/s 
Chirp repetition period Tr 74us 
Number of chirps Ne 240 
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Figure 5.10: Measurement results of the respective one-dimensional (a) range FFT for all N¿ chirps 


overlaid and (b) Doppler FFT for all N, samples overlaid bypassing the fractional delay 
unit O 2021 IEEE [DANZ21] 


response filter with Nyap = 19 taps was implemented on the RTS back end. Its 
adaptive design allowed the filter coefficients to be stored in the DDR memory 
attached to the FPGA and to be loaded at run time. The delay of the virtual 
target was adjusted every Ta = T./2 = 37 us, requiring a delay increment of 


2ARy  2Tavsi 
Ata = u EA 12,33 ps (5.38) 
co co 


This update period was chosen in accordance to the explanations in section 
5.4.2. From this, the number of different fractional delays and thus the number 
of filter coefficient sets that were needed to generate all inter-sample delays can 


be calculated to 
AR co 


N; Z= —— Z _ — — 
AT AR, DB 


zx 162 (5.39) 
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Figure 5.11: Measurement results of the respective one-dimensional (a) range FFT for all N. chirps 
overlaid and (b) Doppler FFT for all N, samples overlaid with the fractional delay unit 
in place © 2021 IEEE [DANZ21] 


In this context, a coefficient set is the collection of filter coefficients required 
for a specific fractional delay. The filter coefficients of each set were calculated 
according to (5.37) with 


Tra = Nfa + ATu , nia € [0, Nfa — 1] (5.40) 


and uploaded to the FPGA's DDR memory. In order to adjust the target delay, 
a new coefficient set is loaded from memory and substituted for the current one 
at the fractional delay filter. This exchange of filter coefficients is clocked at the 
update rate Tą described earlier. The results of the measurement are depicted in 
Fig. 5.11, where the range and Doppler profiles are shown. It can be observed 
that the target migrates of several bins which verifies the successful synthesis 
of the range-Doppler migration phenomenon. The number of visited bins can 
be calculated with (5.9) to Nmig = 3, which stands in good agreement to the 
measurement results when considering the uncertainty of the exact position of 
the initial range and Doppler peak within the occupied bin. 
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5.5 Practical Considerations 


As was demonstrated in this chapter, the synthesis of the migration phenomenon 
was successfully realized and demonstrated on a digital radar target simulation 
system. However, when implementing the proposed approach for a commercial 
solution a few practical considerations have to be taken into account. 


First, the filter coefficient sets that were utilized for the measurement were 
generated for the particular velocity of the virtual target. Different velocities 
require different delay increments, thus different coefficient sets must be gen- 
erated. These must be pre-computed and stored in the system's memory, since 
calculating them in real time is impractical given the required update frequency. 
One possible solution to this is to compute the coefficient sets for the smallest 
detectable velocity and use it as a LUT with an incremental address input whose 
step size can be adjusted according to the velocity of the virtual target and thus 
to the required update frequency. Alternatively, the fractional delay filter can 
be designed as a polynomial filter using the Farrow structure whose tunabil- 
ity allows simple handling of the coefficients and eliminates the necessity to 
recalculate them [Smi10]. 


Next, the hardware resources used to implement the fractional delay filter must 
be considered. The number of filter taps Map correlates with the number 
of complex multipliers and thus with that of DSP slices needed. As can be 
concluded, a compromise between quality and quantity of the fractional delay 
filters, meaning the filter response flatness and number of filters, must be made 
in order to deal with the limited resources of an FPGA. A configurable target 
generation unit design with a demand-driven dynamic allocation of the filters 
only in RTS channels with moving radar targets that show migration behavior 
could tackle this problem. 


Finally, the adaption of the approach towards analog radar target generation 
demands an analog implementation of the required delay increments. Possible 
solutions are true-time-delay phase shifters [TMH*03], surface acoustic wave 
(SAW) filters [Hol74], bulk acoustic wave (BAW) filters and co-fired ceramic 
low-temperature (LTCC) that need to be controllable [ALP*17]. 
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In the past, there have been essentially three techniques to simulate and control 
the angle of arrival of a virtual radar target generated by an RTS system. 
The most intuitive method is to adjust the physical angular position of the 
RTS channel according to the desired AoA of the target. This can be done 
by mechanically moving the respective FE along a circular trajectory that is 
centered on the position of the RuT [GR17, ABB*21,KEY, Way11]. However, 
this comes with the disadvantage that only a limited number of FEs, and 
consequently different angles, can be realized since they block each other’s 
angular movement. Furthermore, the simulatable maximum lateral velocity of 
the targets is limited by the traverse speed of the mechanical setup. 


The second method remedies this deficiency by rotating the RuT itself rather 
than the FEs [dSP22, Kon21], which allows much quicker lateral target move- 
ments. This technique, however, was developed for the application in end- 
of-line radar calibration tests and thus is only capable of simulating a single 
independent target angle and is not suited for integrated ADAS validation test 
procedures. 


Lastly, for the third method, which this work focuses on, multiple FEs are dis- 
tributed in the lateral domain and the control of the AoA of a virtual target is im- 
plemented by electronically switching between them [Key22, Roh22, VGS*21]. 
It enables the utilization of a high number of FEs and is therefore the only 
method capable to represent complex traffic scenarios as they are needed for 
thorough ADAS validation testing that includes current or future radar sensors. 
The approach allows fast lateral target movement and a large amount of virtual 
targets simultaneously, but comes at a high price as it requires many FE units 
whose hardware components are costly. Furthermore, it inevitably introduces 
angular deviations of the simulated AoA from its set-point, as targets can only 
be generated at the discrete physical angular positions at which the FEs are 
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located. Each virtual radar target is synthesized by a single RTS channel and 
its AoA is defined by the physical angular position of the corresponding FE. 
Consequently, the rasterization of the simulatable angle of arrival of the RTS 
system is driven by the lateral spacing of the FEs A@. The maximum angle error 
occurs for a virtual radar target whose AoA ought to reside exactly between 
two neighboring FEs and thus equals half the FE lateral spacing (a. = A@/2). 
Although this error can be reduced by decreasing the angular spacing of the 
FEs, this comes at the expense of an increased number of FEs in order to cover 
the same angular view span, may cause inter-pair coupling between the trans- 
mit and receive antennas of the RTS channels, and is further limited by their 
physical dimensions. 


The new approach that is presented in this chapter overcomes these shortcom- 
ings and enables the simulation of virtual radar targets at an arbitrary angle of 
arrival. It is based on the method of lateral distributed FEs, which will serve as 
a reference for the later evaluation and will be referred to as the conventional 
approach in the following. By superimposing the echo signals of adjacent FEs, 
this new technique allows to drastically reduce the introduced angle error while 
at the same time reducing the number of FEs needed. It is neither limited in 
regards of the methodology of the RTS, as it is applicable for analog and digital 
systems, nor by the modulation scheme of the RuT. 


In the further course of this chapter, first, the basic concept of the new approach 
and its advantages are explained. Next, an extended signal model based on 
that of a generic RTS system in chapter 4.6 is developed that facilitates the 
comprehension of the method and its limitations. Thereafter, the constraints 
of the approach and the calibrations required to satisfy them are elaborated on. 
Subsequently, the results of a measurement campaign that verify the functional 
operation and demonstrate the accomplished linearity of the technique are 
presented. After that, investigations regarding the angle accuracy that can be 
achieved in relation to the lateral spacing of the utilized FEs and the coherence 
properties of the respective RTS channels are exhibited, which allows to draw 
conclusions concerning the number of FEs needed for the overall setup. Next, 
the quasi-monostatic character of RTS systems, due to the physical separation 
of the receive and transmit antennas of the FEs, and its influence on the angle 
accuracy are analyzed. Finally, the approach is extended for the simulation of 
virtual targets in both angular domains simultaneously which makes it suitable 
for the application of 4D radar sensors that are capable to resolve targets in the 
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azimuth and elevation plane concurrently. A further measurement campaign 
verifies the successful implementation of the extension. 


6.1 Superposition of Adjacent Echo Signals 
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Figure 6.1: Concept of the superposition of the returning radar signals of two adjacent RTS channels 
(red and yellow) forming a composed virtual radar target (blue) © 2021 IEEE [DNPZ22] 
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The basic concept of the new approach is to superimpose the returning radar 
signals of two adjacent RTS channels to form a composed virtual radar target 
whose AoA resides between the physical angular positions of the respective 
RTS FEs and that can be adjusted arbitrarily. The artificial target’s AoA can be 
controlled by the amplitude attenuation of the respective RTS channels. The 
concept ofthe approach is visualized in Fig. 6.1. It allows to strongly reduce the 
angle error of an RTS system, as it provides the capability to synthesize virtual 
targets at an arbitrary angle of arrival that lies between the corresponding FEs 
and at the same time decreases the number of FEs needed. In the conventional 
approach, the amount of FEs required could be determined by dividing the RuT’s 
field of view by the desired angle accuracy, or respectively the FE spacing of the 
RTS. In contrast, with the new approach the desired angle accuracy is mostly 
independent of the FE spacing and thus the number of FEs only correlates with 
the RuT’s field of view divided by its angular resolution. Detailed investigations 
regarding the achievable angle accuracy and the required number of FEs can 
be found in chapter 6.6. 


6.2 Extended Signal Model 


For an analytical investigation of the approach and to identify its constraints, 
an extended signal model was developed that is based on that of a generic RTS 
system in chapter 4.6. It demonstrates the capability of the concept to control 
the synthesized AoA through the signal attenuation of the respective RTS 
channels. The model focuses mostly on the signal phase since it is essential for 
the successful superposition of the respective radar signals to form a common 
virtual target, as will be shown later. For the sake of simplicity, a frequency- 
modulated continuous wave (FMCW) radar will be assumed. Nonetheless, the 
underlying principle of this approach operates independently of the modulation 
scheme of the radar. 


Redefining the angular position of an RTS front end 6,, formerly introduced in 
(4.44), to be specific for a particular RTS channel q with q € [1, O], allows the 
phase shift of the antenna element na of the RuT to be expressed as 


da sin(@q) 
A AAAAAAA<<< KO NA 


Pa =2n (6.1) 
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where da denotes the antenna element spacing. Thus, the free space signal 
delay between a particular RTS channel q and a particular antenna element na 
can be described as 


2Re,g A da sin(@q) l 


co co 


Tog = Tix + Trx = na (6.2) 
where R.,g is the physical distance between the specific FE and the phase 
center of the RuT, which are assumed to be equal for the transmit and receive 
paths. Furthermore, the signal attenuation and the artificial delay of a specific 
RTS channel are redefined as A, and Trts,q, respectively. These adjustments 
in relation to the generic RTS signal model allow the radar signal after the 
beamforming process to be expressed as 


Sa,g(@) = Ag © NANN; - sinc (2 [sin(@,) - sin(a)] ma) (6.3) 


, B B da sin(@,) Na — 1 
- exp fiar lr + z) To,q + (Fs + z) Trts,q + er 


with the substitute for the invariant delay 
T0,q = —_ (6.4) 


Comparing this expression to that of the generic model in (4.46) unveils that 
their only difference is that the phase and the shift of the sinc-function are RTS 
channel specific. If the virtually generated radar targets of two neighboring 
RTS channels (Q = 2) are detected in the same range and Doppler bin, their 
complex-valued superposition can be described as 


Q 
Sala) = X sala) = sa (a) + 4,2(2) (6.5) 
q=1 
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which in turn can be written out in full as 


Sala) =NaNeN; A - sinc (> [sin(@;) - sin(a)] ma) (6.6) 


i B B da sin(@;) Na — 1 
"exp {ian (+ 3) To,1 + ES 3) Trts,1 + 2 ( 1) = | 


A 2 


+ A> - sinc (2 [sin(82) - sin(a)] Na) 


: B B da sin(@2) Na -1 
"exp izr (fs ae 3) To2 + Im T 3) Trts,2 + 2 ( 2) = | 


A 2 


This complex-valued superposition of two sinc-functions and their respective 
phases compose a common maximum that lies in between their individual 
maxima. Fig. 6.2 plots this superposition (SP) as an example for two different 
configurations of the RTS channel attenuations A; and Az. As can be seen in 
the diagrams, the maximum of the superimposed function can be controlled by 
the amplitudes of the individual channels and reveals the AoA of the commonly 
formed virtual radar target. In order for the individual signals to superimpose 
additively, their phases must be coherent, which is explained in greater detail in 


Normalized amplitude 
Normalized amplitude 
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Figure 6.2: Simulation of exemplary superpositions (SPs) of two adjacent echo signals forming a 
common virtual radar target for two different amplitude configurations with Na = 8 
and da = A/2 
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chapter 6.3. In the upcoming section, the quantified relationship between the 
individual signal amplitudes and the resulting AoA of the superimposed radar 
target are presented. 


6.2.1 Control of Simulated Angle of Arrival 


Determining the maximum of the commonly formed function SaA(a) for a 
given amplitude configuration is not analytically straightforward, since the 
corresponding variable a is nested in two independent sinc-functions. However, 
to control the simulated angle of arrival only requires to calculate the amplitude 
configuration needed for a given target angle a. For this, Sa (œ) is derived 
according to œ and set equal to zero 


3Sa (a) _ Osa, (q) A Osa, 2 (a) 7 


da da da y en 


Since the sinc-terms of the function are only connected commutatively with 
each other, their derivatives can be formed individually and substituted with 


ô sinc (4 [sin(6,) - sin(@) | Na) 


8q(@) = Fa (6.8) 
cos(q@) sin (<4 (sin(@,) — sin(a)) Na) 
5 na [sin(6,) - sin(a) | Na 
cos(@) cos (4 [sin(6,) - sin(a)] Na) 
E sin(@,) — sin(a) 
allowing the generic individual function’s derivative to be expressed as 
Mralo) = AgNaNeNs ` 8q(@) (6.9) 


: B B 
‚exp fr lr T 3) T0,q + [m + 3) Trts,q 


„Fa sin(@,) Na —1 
A 2 
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Solving (6.7) for the relation of the RTS channel amplitudes 41/42 leads to 
da Na- 1 
A 2 


+ (fs + z) (70,2 - 70,1) 


+ (Fs + 3) (Tres, 2 = 61) ) 


which can be utilized to determine the channel attenuations ratio for a desired 
angle of arrival ag, of a virtual target. To address the remaining underdetermi- 
nation of the individual attenuations, the amplitude-based superposition of the 
signals to a common virtual radar target with a given RCS can be used 


Aj _ _ 82( drat) li ( sin(62) = sin(81)) (6.10) 


Az 7 1 (Gigi) 


(6.11) 
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Figure 6.3: Simulation of the amplitude attenuations of the respective RTS channels over the set- 
point angle with an antenna configuration of the RuT of Na = 8 and da = 1/2 
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where P is the power of the superimposed radar signal and Ra the ohmic 
impedance. Fig. 6.3 shows the amplitude attenations of the respective RTS 
channels as a function of the set-point angle aigt. 


6.3 Constraints 


The successful constructive interference of the individual signals is subject to 
certain constraints. As mentioned before, the individual radar targets must 
be detected in the same range-Doppler bin, in order for the signals of the 
spatial-spectrum to be superimposed in the subsequent beamforming process. 
Additionally, the spacing of the FEs A@ = 62 — | and the corresponding dis- 
tance between the individual signal maxima after beamforming must be small 
enough so that the commonly formed signal features only a single maxima. 
Otherwise, two individual maxima emerge which may be unintentionally de- 
tected as two angularly resolvable targets. This case is depicted exemplarily 
in Fig. 6.4. As can be concluded, the threshold of the FE spacing for a suc- 
cessful superposition correlates with the width of the sinc-function, and thus 
with the RuT’s angle resolution. Further research on this topic is presented in 
chapter 6.6. 
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Figure 6.4: Simulated unsuccessful superposition (SP) of two adjacent echo signals forming two 
angular resolvable radar targets due to excessive FE spacing with Na = 8 and da = 4/2 
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The next constraint is that two adjacent RTS channels can only simulate a single 
common target per range-Doppler bin, since the signals of multiple individual 
targets would be coherent in regards of phase and frequency and therefore 
overlay additively and form a composite target at an intermediate AoA. 


Phase Coherency 


A major constraint concerns the individual signals’ phases. For a purely con- 
structive superposition, the phases of the signals must be coherent, or else the 
signals will partially or completely cancel each other out, resulting in a distorted 
angle detection. For this, the individual phase controlling terms in (6.6) must 
be set equal leading to 


exp {iva.i} = exp {iva} (6.12) 
with the generic phase 


da sin(6,) Na — 1 
A 2 


B B 
Pag = 20 lr + 2) T0,q + (Fs + 3) Trts,q + (6.13) 
From (6.12), it can be concluded that any phase difference between the two 
individual signals must be an integer multiple of 27 


Ag = Pa - Yan = nn, neN (6.14) 


The individual signal phases are very sensitive to variations in the physical 
distance R¿ y between the RuT and the respective FE, hence even small inac- 
curacies in the mechanical mounting of the FEs AR. = R¿ 1 — R.,2 lead to a 
significant relative phase offset and potentially destructive interference of the 
individual signals. For example, assuming a carrier frequency of fer = 77 GHz, 
a relative radial position offset of one of the FEs of only AR. = 4/4 ~ 967 um, 
causes an overall phase shift of Ay = 7 and the extinction of the composite 
signal. This is shown exemplarily in Fig. 6.5, where the detected target angle is 
plotted as a function of the set-point angle for different radial distance deviations 
AR, and phase shifts Ay. As can be seen, a distance offset of AR. = 1/2 does 
not affect the detected AoA, since it corresponds to a phase shift of Ay = 27, 
which in turn has no influence on the phase coherency. Therefore, it appears as 
the radial distance of the FEs must only match in multiples of 2/2. 
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Figure 6.5: Simulation of the detected angle of arrival for different (a) radial distance offsets AR. 
and (b) phase shifts Ay with Na = 8 and da = 4/2 


However, if the distance discrepancy AR. is larger than just fractions of the 
wavelength, the positions of the individual maxima within the FFT bin in the 
range profile deviate from another. This causes the individual amplitudes to 
differ at the center ofthe FFT bin, which is used for the subsequent beamforming, 
leading to an error in the detected AoA, as can be seen in Fig. 6.6. Since the 
distance offset AR. in this plot is chosen in multiples of 5A = 19 mm, this does 
not result in any phase shifts. The occurring angular error is similar to that 
caused by mere amplitude deviations AA = Aı — Aa. Since the AoA of the 
commonly formed radar target is controlled by the amplitude ratio of the RTS 
channels, any deviations in their amplitude levels will result in errors in the 
detected AoA. Therefore, the signal amplitudes of the involved channels must 
be balanced. For illustration, Fig. 6.6 also shows exemplarily the detected over 
the set angle for different amplitude deviations AA. 


Since the RTS channels are under the influence of various hardware imperfec- 
tions and positioning inaccuracies, it becomes clear that in order to meet these 
constraints and thus reduce the occurring angle error, the channels must be 
calibrated both in conjunction and individually with regard to delay, phase and 
amplitude. 
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Figure 6.6: Simulation of the detected angle of arrival for different (a) radial distance offsets AR. 
and (b) amplitude deviations AA with Na = 8 and da = 1/2 


6.4 Calibration 


As pointed out in the previous section, differences between the RTS channels 
in terms of delay, phase and amplitude cause errors in the angle detection and 
thus must be calibrated out. A prerequisite of the RTS system as a whole, 
and thus also for the calibration, is that it must be applicable for commercial 
radar sensors, meaning that only the data of the detected point targets after 
the range-Doppler FFT, the beamforming, and the subsequent detection is 
available [Con21,Rob21]. That is the target range, velocity, RCS and AoA. 


The calibration procedure is executed in three steps, which are shown in a 
flow chart in Fig. 6.7. First, the physical angular positions of the FEs are 
measured individually. Next, in order to meet the criterion that the individual 


Measure FE Range & Amplitude Phase Calibration 
angle —>| Calibration I+ sweep: Ay 
measure: 0, adjust & monitor: Tg & A, monitor: angle error 


Figure 6.7: Calibration procedure flow chart 
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radar targets must be detected in the same range-Doppler bin, the target ranges 
of the respective RTS channels are measured individually and calibrated by 
adjusting the artificial delay 7;ts,¿. The same can be performed for the amplitude 
respectively RCS calibration using the channel attenuation A,. Finally, a 
collaborative phase calibration is executed. 


Phase Calibration 


The calibration of the signal phase and the delay within an FFT bin presents 
itself as not as simple, since they are not output directly by commercial sensors. 
Furthermore, even if they were, such a calibration would require the signal delay 
respectively target range to be determined with an extremely high accuracy. 
However, high-precision radar range estimation down to the order of fractions 
of the wavelength A, as it is required here, demands for a high signal-to-noise 
ratio (SNR) and can nevertheless only provide a relative rough estimate of the 
actual delay together with a statistical specification for its estimation accuracy, 
as has been shown in [Sch17a]. Therefore, the calibration is realized on the 
basis of the angle error, whose sensitivity to delay and phase differences can 
be used. Even though this method bears a certain indeterminacy and does not 
allow for an exact determination of the signal delay and phase, it operates in 
the sense of the goal of the calibration to minimize the relative phase offset and 
thus the angular error of the superimposed target. 


The calibration is performed by monitoring the occurring angle error while 
sweeping the signal phase of one of the RTS channels and keeping that of 
the other constant. The angle error is calculated as the absolute value of 
the difference between set-point and actual AoA. The set-point can be chosen 
freely between the two respective FE positions and the required RTS channel 
attenuations are determined according to (6.10). The phase sweep can be 
executed using a phase shifter or, as (6.13) reveals, using the artificial target 
delay Tris, whose adjustment causes a phase shift corresponding to 


OPa,q 
OTrts,q 


=2n (m 4 2) (6.15) 


For illustrative purposes, Fig. 6.8 shows exemplarily the simulation results of 
the calibration process plotting the calculated angle error as a function of the 


105 


6 High-Precision Angle of Arrival Simulation 


5 T bu T I T 5 T T T T T 
— Phase shift AYrts — Phase shift Ayyts 
4 —— Delay shift Ants A = Delay shift Arts 
20 50 
Q © 
I =) 
23 2 3 A 
E E 
= 2 
3 3 
2 2 2 2 + 
eb ab 
g =] 
< < 
1 1F 
0 | | | | | 0 | | | | | | 
—2-15-1-050 05 1 15 2 —2-15-1-050 05 1 15 2 
Phase offset in 7 Phase offset in 7 
(a) AR, = 4/8 (b) AR. = (104 1/8) A 


Figure 6.8: Simulation of the phase calibration with the angle error over the phase offset using 
phase shifts Apis and delay shifts At,, for two different physical range offsets AR. 


phase shifts induced by a phase shifter Apis € [-27, 27] and a delay offset 
ATrts € [—1 ns, Ins] for two different initial physical distance offsets AR¿. As 
can be seen, in both plots the angle error arrives at a minimum for a phase shift 
of Ay = -n/2, which is in accordance with the previous elaborations since it 
eliminates the initial phase shift introduced by AR¿. However, as is noticeable 
in the comparison of the two graphs, there are differences in the minimum 
values of the angle error for the respective parameter sweeps. This is due to the 
difference between the two initial distance offsets of 104, as it does not cause 
any phase shift, but it does cause a deterioration of the amplitude balance which 
in turn introduces an angle error as mentioned before. The difference of the 
minimum angle error between the two different phase sweeping methods within 
one graph can be explained by the range shift caused by the delay shift Arııs, 
which also leads to an amplitude imbalance, but of much smaller magnitude. 


It is for this reason, that neither of the two phase calibration methods clearly 
outperforms the other, given the uncertainty of the exact initial position of the 
individual target peaks within the range bin, as it is to be expected in a real 
measurement environment. In order to nevertheless make a selection between 
the two methods, their implementation complexity in the respective system 
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topology can be considered. Either method is fairly easy to realize in the 
digital domain with complex multipliers. On the other hand, the implemen- 
tation of a variable, frequency-invariant phase shift for an analog signal with 
a high relative bandwidth is rather difficult when only analog components are 
available. The analog implementation of a small-step delay shift, while also 
non-trivial, is more viable [TMH*03,H0174, ALP*17]. Therefore, the following 
verification measurement will focus on the true-time-delay sweep calibration 
method in order to demonstrate not only the practical functionality of the signal 
superposition approach, but also the feasibility of its implementation. 


6.5 Verification Measurement 


In order to verify the functioning of the proposed approach, a measurement 
campaign was conducted and published in [DNPZ22]. For this, the RTS 
system that is described in chapter 4 was set up in front of a Texas Instruments 


BB RTS front ends 


Figure 6.9: Photograph of the arbitrary angle of arrival measurement setup O 2021 IEEE [DNPZ22] 
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Table 6.1: Radar parameters 


Name Symbol Value 


Carrier frequency fa 77 GHz 
Bandwidth B 1 GHz 
Chirp repetition period Te 41.33 us 
Number of chirps Ne 120 
Sampling rate fs  25MS/s 
Number of transmit antennas N 2 
Number of receive antennas Nix 4 
Antenna element spacing da A/2 


AWR1843BOOST radar sensor, that was utilized as the RuT, at a physical 
distance of R. = Im. A photograph of the setup of the measurement is shown 
in Fig. 6.9. The receive and transmit front ends of the RTS were placed in a 
semi-circle formation, vertically aligned and behind a metal sheet that was used 
to block off undesired radar reflections of any structures in the background. 
The carrier frequency of the RTS was set to fits = 500 MHz and the sampling 
frequency of the ADCs and DACs to forts =4GHz. The radar sensor was 
configured according to the parameters in Table 6.1. 


Calibration Measurement 


As a first step, a calibration measurement in accordance with the method de- 
scribed in the previous chapter was conducted. For this purpose, two RTS 
channels were synchronized in regard of their amplitude and the range-Doppler 
bin occupied by the virtual radar target that they generate. For this individual 
measurements were performed with only one channel active at a time. The 
front ends of the channels were solely distributed in azimuth, with no elevation 
differences and positioned at the azimuth angles of 0; = 3.4° and 62 = 12.2° 
from the perspective of the RuT. The phase calibration was performed by moni- 
toring the angle error, which is calculated as the difference between the detected 
and set-point angle, keeping the delay of the first RTS channel 7;:s,ı constant, 
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Figure 6.10: Simulation and measurement for the phase calibration with the angle error over the 
artificial delay offset ATrts © 2021 IEEE [DNPZ22] 


while sweeping that of the second channel T;1s,2 = Trts,ı + ATs in an interval of 
ATzıs = [-0.5 ns, Ins] with a step size of 25 ps. The smallest delay increment 
realizable with sampling buffering and the sample frequency used to discretize 
the signal ATyts bur = 1/ forts = 0.25 ns exceeds the delay step size required for 
the phase calibration, which is why the fractional delay filter unit that is de- 
scribed in 5.4.3 was used. Fig. 6.10 shows the measurement of the calculated 
angle error as a function of the artificial delay offset Ar,ıs. As can be observed, 
the angle error arrives at a minimum for AT;ıs = 0.1 ns, which is used to calibrate 
the RTS channels for the subsequent measurements. In addition, the simulation 
of the angle error is plotted in order to verify the analytical derivations of the 
extended signal model in section 6.2. Measurement and simulation show good 
agreement, as only small deviations occur that can be traced back to differences 
between the simulated and the actual physical positions of the front ends in re- 
gards of range and angle. For the simulation, the positions were determined by 
a separate measurement, whose inaccuracies lead to the mentioned differences 
in simulated and actual FE positions. 
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Linearity Measurement 


Next, the achievable linearity and its corresponding error of the AoA of the 
superimposed virtual radar target were measured. For this purpose, the RTS 
channels were synchronized in regard of amplitude, delay and phase according 
to the previously performed calibration measurement. For the linearity mea- 
surement itself, the set-point angle was increased in 100 steps in the interval 
of aigt € [01,02] = [3.4°, 12.29] and the required channel attenuations were 
calculated according to (6.10) and (6.11). In Fig. 6.11, the detected AoA over 
the set-point angle ag is depicted. Additionally, the results of a simulation 
of the detected angle according to the signal model in section 6.2 are plotted 
which show good agreement with the measurement results. In Fig. 6.12, the 
calculated angle error of the measurement and the simulation can be seen. Al- 
though the angle error of the measurement exceeds that of the simulation, the 
maximum measured error of less than 0.18° correlates to only 2.1% of the 
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Figure 6.11: Simulation and measurement of the detected angle over the set-point angle O 2021 
IEEE [DNPZ22] 
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Figure 6.12: Simulation and measurement of the angle error over the set-point angle © 2021 
IEEE [DNPZ22] 


angular separation of the RTS FEs (A@ = 8.8°). This represents a reduction of 
the angle error by 


0.18° 


= 205: 1 
50% Ad 95.9 fo (6.16) 


when compared to the conventional approach, in which no signal superposition 
1s performed and virtual targets are only simulated by a single RTS channel each. 
Furthermore, it can be assumed that the angle error occurs due to a remaining 
amplitude and phase offset between the utilized RTS channels due to a non- 
ideal calibration. The latter could be further reduced by an iterative calibration 
with decreasing step sizes for A7;ts, but can not be eliminated completely as 
the aforementioned restrictions still apply. The simulation again serves as a 
reference for the underlying theory developed in the previous chapters. 


Dynamic Multi-Target Measurement 


For the final verification measurement, a scenario with four dynamic virtual 
targets was generated by the RTS. All targets were created by the same two RTS 
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Table 6.2: Virtual Target Properties 


Target Range Velocity Angle 


1 335m  0m/s 7° 
2 37m 4m/s 4° 
3 45m -2m/s 10° 
4 52m -—5m/s 11° 


channels simultaneously, but with different range, radial velocity and angle of 
arrival characteristics which are stated in Table 6.2. For this, the channels were 
calibrated in the same manner as before and also the Doppler shift components in 
the target generation unit were synchronized in regard of their phases. Fig. 6.13 
illustrates the range-Doppler and range-angle plot with red crosses indicating 
the detected targets of the measurement. It can be observed, that the targets 
are generated with the correct range, velocity and angle of arrival features. 
This proves the functionality of the approach to generate dynamic, multi-target 
scenarios, which are required for the validation testing of ADAS functions and 
the radar sensors involved. The spurious peaks that occur in both plots with 
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Figure 6.13: Measurement of a dynamic multi-target scenario with red crosses indicating the de- 
tected targets O 2021 IEEE [DNPZ22] 
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an additional range offset to the intended targets are caused by the mismatch 
of the characteristic impedance between the front ends and the coaxial cables 
connected to them. The radar signal travels back and forth between the back 
and front end of the RTS, creating ghost targets with range offsets of multiples 
of the cable length. The static reflection in close vicinity to the radar can be 
assigned to the mechanical structure setup of the RTS. 


6.6 Achievable Angle Accuracy 


In the previous chapter, the practical functionality of the presented approach 
for generating superimposed virtual radar targets at an arbitrary angle of arrival 
and thus its capability to increase the angular accuracy of an RTS system 
was demonstrated. However, the verification measurement was conducted for 
only one specific front end separation value. In order to further investigate the 
ability and limitations of the approach to reduce the angular error of a simulated 
target and thereby to estimate the number of RTS channels required, this section 
examines the achievable angle accuracy as a function of the front end separation 
A8 and the respective channel imbalances. In this context, angle accuracy can 
be equated with the angle error with which the virtual targets are generated 
by the RTS. This allows an in-depth analysis of the aforementioned constraint 
that the front end spacing must not exceed a certain threshold in order for the 
RTS channels to form a common virtual radar target. The elaborated findings 
provide a basis for the systematic decision making regarding the total number 
and placement of the RTS FEs. 


For this purpose, another series of measurements was performed using the 
same RuT and RTS system with the same configurations as in the verification 
measurement. Two front end receive-transmit pairs were positioned at a distance 
of R¿ ~ 1m to the RuT. The first pair was mounted statically on a tripod, while 
the second one was installed on the end of a wooden beam that was attached to 
a rotation motor that resided below the RuT. This setup allowed to change the 
angular position of the second FE and thus to adjust the FE separation angle 
A98 at will. A photograph of the assembly is shown in Fig. 6.14. 


For the actual measurement process, the FE separation angle A@ was increased 
stepwise. In each step, the RTS channels were calibrated with regard to range, 
amplitude and phase, and subsequently a superimposed virtual radar target 
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RTS front ends 
= 


Figure 6.14: Photograph of the achievable angle accuracy measurement setup 


was created whose intended AoA as: was swept from the position of the first 
FE to that of the other in Nsp = 100 steps, in the same way as in the linear 
measurement procedure. For each step the relative angle error was computed as 
the absolute value of the difference between the respective set-point @tgt(Mstep) 
and detected angle @det(Nstep) over the front end separation angle A0 


Qtgt (Nstep) = det (Nstep) 


m (6.17) 


Qe,rel (Nstep) = 
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where "step € [1, Nstep] describes the step index. Then, for each FE separation 
step the maximum relative angle error was calculated among all simulated angle 
steps 


Qe rel, max (AQ) = max We rel (step) (6.18) 


Mstep 
Next, the front end separation AQ was increased and the procedure was repeated. 


The results of the measurement together with those of a simulation that is 
based on the extended signal model developed in section 6.2 are shown in 
Fig. 6.15, where the maximum relative angle error is drawn as a function of the 
FE separation A@. This form of presentation was chosen to liberate the results 
from the parametric configuration of the RuT (angle resolution) and to provide 
a generic representation of the capability of the approach. The simulation takes 
into account the occurrence of amplitude (AA) and phase mismatches (Ay) 
between the RTS channels, which can be caused by an imperfect calibration or 
potential hardware influences. These mismatches are potentially also present 
in the measurement results, but to an unknown extent for which the simulation 
serves as a reference. Due to the geometrical dimensions of the FEs the 
minimum FE separation realizable for the measurement was 4°. A horizontal 
dashed line indicates the maximum relative angle error of the conventional 
angle simulating approach for RTS systems with lateral distributed FEs (50 %) 
as a reference. 


The results show that the angular error of the system becomes more sensitive 
to amplitude and phase imbalances between the RTS channels for larger FE 
separations. In particular, the simulated phase deviation of Ay = 7/2 causes a 
partially destructive interference of the superimposed target and the occurence 
of an unintended second target, thus a rather strong distortion of the detected an- 
gle. This also complies with the previous simulation results shown in Fig. 6.5a. 
The results underline the importance of the preceding calibration and demon- 
strate the dependency of the angular accuracy thereon. In general, measurement 
and simulation are in good agreement, which verifies the signal model and its 
ability to estimate the angle error to be expected given the separation of the FEs 
and the imbalances of the RTS channels. Conversely, the model can also be 
used to determine the number of FEs needed when the field of view (FoV) of 
the RuT, the RTS channel imbalances and the angle error tolerable are given. 
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Figure 6.15: Measurement and simulation of the achievable angle accuracy as a function of the front 
end separation A@ in the presence of (a, b) amplitude and (c, d) phase imbalances 
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Both estimations can support the design and layout process of high-precision, 
angle simulating RTS systems. 


Angular Boundary Condition for Superposition of Echo Signals 


The maximum relative angle error rises sharply for a FE separation that is too 
large for the superimposed signal to form a single, common peak. For this FE 
spacing and onward, two distinct peaks are formed in the angular profile and 
the AoA estimated by the radar deviates significantly from its set-point. Coun- 
terintuitively, this turning point does not coincide with the definitions of the 
angular resolution (half-power beamwidth Aa, and Rayleigh criterion Aa...) of 
the utilized radar sensor given in section 2.5.1 and marked with dashed, vertical 
lines in Fig. 6.15. It lies beyond either of them, as can best be witnessed in the 
case of the ideal simulation (AA = 0 and Ay = 0), where the maximum relative 
angle error increases sharply only for a FE separation not negligibly larger than 
the two. This discrepancy is a result of the different coherence assumptions 
between the superposition of radar echoes with the intent of composing a virtual 
radar target at a controllable AoA and the ability of a radar to discriminate two 
independent targets in the angular domain, that is, the angular resolution. 


Both definitions of the angular resolution presume two spatially separated tar- 
gets of equal range, radial velocity and electromagnetic size whose reflected 
radar signals are incoherent in regard of their phases. Due to RCS fluctuations, 
the individual echo signals exhibit a random phase relation, which can lead 
to constructive or destructive inference of the individual targets within a sin- 
gle measurement. Throughout the course of multiple measurement cycles, the 
relative phase of the individual signal is uniformly distributed, hence the super- 
position of the individual target signals according to the angular resolution is 
considered as an incoherent summation only in terms of the signal power, that 
is, the absolute value of the amplitudes and neglecting any phase information. 
Therefore, the angular resolution can be thought of as a statistical expression 
for the separability of two targets in angular domain, that is not applicable to 
the superposition of two phase-coherent echoes as encountered in this work. 


In contrast, the concept of superimposing the returning signals from two ad- 
jacent RTS FEs for the purpose of generating a virtual radar target with an 
arbitrarily controllable AoA as presented here requires not only phase co- 
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herency, i.e., a constant phase relation of the RTS channels involved, but also 
a phase offset of multiples of 27 (Ay = pı - ya = 2nr ,n € N). This con- 
straint is explained in more detail in section 6.3 and the calibration process 
discussed in section 6.4 is employed to satisfy it. The result is a phase coherent, 
complex-valued superposition of the individual radar echoes with respect to the 
individual signal amplitudes and phases. 


The applicable boundary condition for the coherent superposition of two radar 
targets to form a common peak, i.e., the coherent angular resolution, is given 
by the first inflection point of the angle profile of an individual target echo. It 
can be found by setting the second derivative of the signal function from (6.3) 
equal to zero 


sa g(a) 
—=0 6.19 
da? ( ) 


and solving numerically for œ, resulting in a coherent angular resolution of 


A 


A = 1.32 
Och Nada 


rad (6.20) 


If the two individual target functions are superimposed at their first inflection 
points, they form a plateau between them, since the function value of the 
first signal decreases from this point onward as quickly as that of the second 
function increases. For a larger separation, two distinct peaks are formed. The 
boundary condition of the coherent superposition is also indicated in Fig. 6.15 
and coincides with the turning point of the maximum relative angle error in the 
case of the ideal simulation. 


The distinction between the coherent and incoherent superposition of two radar 
echoes is illustrated exemplarily in Fig. 6.16 for the borderline case of a angular 
separation of Aare < A0 = 18° < Age. As can be seen, for this particular 
separation angle the incoherently superimposed signal already features two 
distinct peaks, whereas the coherent superposition still forms a single common 
one. The first inflection points of the individual signals are marked by dots. 


Therefore, it can be concluded that it is necessary to take into account the phases 
of the RTS channels whether the here presented concept to superimpose virtual 
targets or the conventional AoA simulation approach is used, since unfortunate 
phase relations among neighboring channels can lead to destructive interference 
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Figure 6.16: Simulated angle profile for coherent and incoherent signal superposition of two radar 
targets after beamforming for a separation angle of A@ = 18° 


canceling out the intended targets. This means that even if the presented 
approach is omitted, a phase calibration or alternatively an RCS fluctuation 
simulation model is required. However, as has been shown, the presented 
concept enables to greatly reduce the number of FEs needed and therefore cuts 
the corresponding hardware costs substantially. 


6.7 Quasi-Monostatic Antenna Displacement 


Another source of angular error in radar target simulation, in addition to the 
previously mentioned FE separation and amplitude and phase imbalances, is the 
quasi-monostatic characteristic of the RTS FEs. The antennas of the RTS for 
reception and re-transmission of the radar signal are placed closely together, but 
still are spatially separated, thus the system can be considered quasi-monostatic. 
Since the simulated AoA of a virtual radar target is, among others, defined by 
the physical position of the respective RTS FEs, this characteristic induces an 
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uncertainty between intended and actual target angle. The radar signal emitted 
by the RuT is received and re-transmitted at slightly different angular positions 
which causes the detected AoA to deviate from the physical position of the 
receiver and transmitter. Fig. 6.17 visualizes the individual signal paths for 
the receive and transmit signal. In radar validation setups where only a single 
angular domain is of interest, this issue can be circumvented by aligning the 
receive and transmit antenna of an FE pair perpendicular to the resolvable 
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Figure 6.17: Conceptual representation of the quasi-monostatic characteristic of radar target sim- 
ulation that leads to differences in signal path lengths and propagation delays. The 
scaling and spacing of the figure were chosen for illustrative purposes and are not 
representative. © 2022 IEEE [DNZ22] 
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angular domain. However, since current and future radar sensors are capable 
of two-dimensional angle of arrival estimation [Con21, Rob21], this technique 
does no longer present a workaround and the uncertainty in the angle simulation 
should be investigated. 


Extended Signal Model 


For this purpose, the RTS signal model presented in 4.6 is adapted and the 
simplifications of equating the angular positions of the receive and transmit 
antenna ofthe RTS FE are removed. The physical position of the FE is split up 
into a receive angle @,, and a transmit angle 6,,. Furthermore, the unification of 
the transmit and receive antenna elements of the RuT to form a virtual antenna 
array is lifted, in order to break down the effects of delay differences in the 
outward and return propagation of the radar signal. For a better illustration, 
the following adapted signal model only considers the spatial separation of 
the receive and transmit antennas in the azimuth domain, as this allows for a 
simplified representation of the model, but still enables the investigation of the 
quasi-monostatic characteristics. With these considerations, the delay of the 
signal due to its propagation through free space to and from the RTS can be 
described as 


Re + dix sin(0;x) 


Tx = E Mix (6.21) 
0 
R. + dix sin(@ 
Tx = Ste -nx (6.22) 


where dix rx specifies the element spacing of the transmit and receive antennas 
of the RuT, respectively, and nıxyrx € [0, Nix/rx — 1] indexes the individual 
antenna elements. The round-trip signal delay can be expressed as 

dix sin (O;x) * Mix + drx Sin (Oix) + Mex 


Te = Tix + Trx = To +t (6.23) 
co 
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with the antenna element independent delay To = 2R./co. With this, the phase 
of the radar signal after range-Doppler processing and detection by the RuT in 
accordance with (4.45) can be formulated as 


Yon (Nix; Nyx) = 20 


fa 3) To + Im z) Trts (6.24) 


y dix Sin(Orx) - Mix + dyx sin (6x) - | 
A 


For the subsequent beamforming, the beamsteering vector must be adapted for 
the non-unified antenna array 


sin(@) 
a 


bip (Mx, Nıx,@) = exp (in (dix + Mix + dıx © nwo} (6.25) 


The spatial spectrum is generated by Fourier beamforming, for which the 
adapted steering vector is applied to the radar signal and the sum is formed 


over all transmit and receive antenna elements and for all angles œ to be con- 
sidered 


sala) = >, > SD (Nix, Mrx) © Dip (Nix, Mx, œ) (6.26) 


Nix-1 Nrx-1 
Mx=- z Mrx=> x 


which can be reformulated to 


sala) = NN, - exp (an 


BY 2R. B 
Ir + z) : E + Im + z) + Trts (6.27) 


€ 


sin(@ 
slas (Mix > 1) + dex (Nix E 1) ) nt ) | 
Nx-1 d 
y exp {ian (sin(0,x) — sin (@)) - ns} 
hx =0 
Nyx-1 


= exp far (sin(x) — sin (@)) “ne 


Nyx=0 
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As before, the partial sum of a geometric series allows to simplify this expression 
to 


: BY 2R B 
sala) = NoN; + exp 4j27 || fe + <= | — + | fits + — | + Tits (6.28) 
2 co 2 
+ dix (Nix 1) sin(6;x) + di (Nix => 1) sin(01x) | 
2A 


- Nix : sinc E (sin(9,x) — sin (@)) Na] 
A dıx . . 
- Nix + sine (= (sin(0íx) — sin (@)) Na] 


This reveals a function that is composed of two sinc-functions, much like that 
in the extended signal model for the superposition of adjacent echo signals in 
chapter 6.2, only this time not as their sum, but as their product. Additionally, 
the widths of the sinc-functions differ from one another and are determined by 
the geometrical dimensions of the receive (dyx Nyx) and transmit antenna array 
(dix Nix) of the RuT. This allows to conclude that the position of the maximum 
of the composite function sa(a) and thus the detected AoA will be located 
between the angular position of the RTS receiver 9,, and transmitter 6x and, 
moreover, depends on the antenna array configuration of the particular RuT that 
is employed. Due to the nesting of œ within the two separate sinc-functions, 
the actual value of the detected AoA can not be derived analytically and must 
be found numerically. 


Measurement 


In order to verify the adapted signal model, a verification measurement was 
conducted using the same RuT and RTS system with the same configurations 
as in the measurements in chapter 6.5 and 6.6. The RuT was positioned at a 
distance of R. x Im in front of an RTS receiver and transmitter. A diagram 
showing the concept ofthe measurement setup can be seen in Fig. 6.18. Multiple 
measurements were performed in between which the spatial separation of the 
RTS antennas was incrementally increased. For this purpose, a linear motor 
moved the RTS transmitter laterally to the RuT while the position of the receiver 
was fixed. An additional rotation motor, that was mounted on the linear motor, 
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A motor 
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Figure 6.18: Concept of the quasi-monostatic antenna displacement measurement setup O 2022 
TEEE [DNZ22] 


compensated for the change in the relative angle of the RTS transmitter as seen 
by the RuT. This was done in order to neutralize any effect that the directional 
antenna pattern of the RTS transmitter might have on the outcome of the signal 
processing chain. The resulting additional range offset was accounted for in 
the simulation, but remained smaller than the range resolution of the radar. In 
Fig. 6.19 a photograph of the measurement setup can be seen. 


The radar signal received by the RuT was processed to estimate the range and 
angle of arrival of the virtual target. For this, a discrete Fourier transform 
followed by a range detection was applied to the signal. Subsequently, the 
AoA was estimated by applying beamforming to the resulting radar signal. The 
results are shown in Fig. 6.20, where the deviation of the estimated AoA from 
the angular position of the RTS receiver is depicted as a function of the lateral 
distance d,ts between the RTS receiver and transmitter. The results of the same 
measurement series were processed with different antenna selections of the RuT 
(2x4, 2x2, 1x4) and plotted in differently colored, thick lines to demonstrate the 
influence of the antenna geometry. The fluctuations in the measurement curves 
are due to multi-path reflections caused by the mechanical structure of the setup, 
as they occurred identically in successive measurements and could be reduced 
by placing absorbers. The output of the simulation is drawn in black and shows 
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Figure 6.19: Photograph of the quasi-monostatic antenna displacement measurement setup © 2022 
IEEE [DNZ22] 


good agreement with the measurement, which verifies the adapted signal model 
and its capability to predict the detected AoA taking into account the quasi- 
monostatic characteristic of an RTS system. Depending on the selected antenna 
configuration of the RuT (Nix X N;x), the signal path differences evaluated to 
determine the AoA of the target are more affected by the number of outward 
(Nix) or return paths (NV, x) of the signal to or from the RTS and are thus more 
dependent on the position of the RTS receiver or transmitter. 


The developed signal model is intended to support the setup and antenna ar- 
rangement process of high-precision angle simulating RTS systems. As has 
been shown, the antenna geometry of the RuT has an impact on the estimated 
angle of arrival, and thus the antenna aperture and orientation of the RuT must 
be known a priori. Usually, they can either be concluded given the parameters 
of the field of view (FoV) and the angular resolution of the RuT, or measured 
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Figure 6.20: Simulation and measurement of the detected angle over the set-point angle © 2022 
IEEE [DNZ22] 


when removing the packaging of the sensor. With this knowledge and the signal 
model presented here, the angle deviation that arises from the quasi-monostatic 
characteristic of the RTS system can be compensated by a minor adjustment of 
the respective FE positions. 


6.8 Two-Dimensional Angle Simulation 


Modern 4D imaging radars are capable to detect targets concurrently in both 
angular domains, azimuth and elevation. Therefore, RTS systems that are 
employed to validate these sensors must be capable to generate virtual radar 
targets at freely controllable lateral positions. To achieve this, the previously 
elaborated approach to superimpose adjacent RTS echoes to form a common 
target at an arbitrary azimuthal position is extended to also cover the elevation 
domain. Through the superposition of the return signal of four neighboring 
RTS channels that are arranged in a vertical square formation, this enhancement 
enables virtual targets to be created at any two-dimensional angular position. 
The arrangement concept of the respective RTS FEs and their superimposing 
radar echoes are shown in Fig. 6.21. 
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In the following, the RTS signal model developed in 4.6 is again extended, this 
time to account for the two-dimensional distribution of both the RuT antennas 
and the RTS FEs. This allows to derive the channel attenuations required 
to control the simulated AoA of the created virtual target in the azimuth and 
elevation domain simultaneously. Subsequently, a measurement campaign that 
verifies the functionality of extension is presented for which, first, the calibration 
process explained in 6.4 is executed for all RTS channels involved, before 
the actual measurement series is performed. Virtual targets at various lateral 
positions distributed in both angular domains are simulated and the detected 
AoAs are compared to their set-points. 


y 


Figure 6.21: Concept of the two-dimensional arbitrary angle of arrival O 2022 IEEE [DNGZ22] 
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Extended Signal Model 


For the extended signal, it is assumed that the transmit and receive antennas of 
the RuT are distributed in only one angular domain respectively, as is typically 
the case in current radar sensors. The Nx transmit antennas are allocated only 
along the vertical axis with an element spacing of dix. The same applies to the 
N,x receive antennas in the horizontal domain with a spacing of d;x. Thus, the 
free space signal propagation delays from the RuT to the RTS FE and back can 
be expressed as 


Reg + dix sin -n 

Tix.g = — = co (a) = (6.29) 
Reg + dıx sin(0¿) cos . 

Trx,q — = l 4) Wa) = (6.30) 


co 


where Yq, 0q € [-90°, 90°] denote the elevation and azimuth angle at which 
the RTS FE is located as seen by the RuT. 


With this, the signal phase after the range-Doppler processing by the RuT 
follows as 


B B 
PD,q (Mix, Mrx) = 27 Ir + 3) To + Im + 5) Tits (6.31) 
„In sin(9¿) cos(Wq) ` yx + dix Sin (Y y) - = 
A 


The two-dimensional steering vector 
271 . ; 
bop (Mx, Mrx, A, B) = exp IT [dıx sin(@) cos(ß) : nx + dix sin(B) : nx] 
(6.32) 


has been derived before in chapter 2.4.2 and is stated again here for the sake 
of clarity. Then the steering vector is applied and the sum is formed over all 
transmit respectively receive antenna elements. Simplifying in the same manner 
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as before using the partial sum of a geometric series, as well as sin(x) = x and 
cos(x) ~ 1 for |x| « 1 leads to 


Nix! Nrx-1 
am I 


sagloB)= Y, DY, Daliwn): bap Mees nap) (6.33) 


Nasi __ Nex-l 
n=- N= 


= AgNeNs - exp {iva.q} 


- Nu  sinc E [sin(y¿) - sin(ß)] Na] 


- Nyx + sinc (= [sin(@,) _ sin(a)| Na] 


with a signal phase of 


Pa q = 20 


B B 
Ir + 3) To + ES + 3) Trts (6.34) 


FA drx (Nix = 1) sin(@q) T dix (Nex = 1) wa 
24 


This time, the sinc-function stretches along both spatial axis and has its main 
maximum at sin(@) = sin(@,) and sin(ß) = sin(,). Through the superpo- 
sition of O = 4 artificial radar echoes and controlling the attenuation of the 
corresponding RTS channels, a common virtual radar target is created whose 
AoA can be steered in the azimuth and elevation domain. The superimposed 
signal after the two-dimensional beamforming can be expressed as 


Q 
Sala, B) = Y saqla, p) (6.35) 


q=1 


sa,1(@, B) + Sa2(a@, B) + 8a,3(@, p) + sarsla, p) 


In contrast to before, deriving the expression above according to œ and £ does 
not allow to determine the required channel attenuations A,, but results in an 
underdefined system of equations due to the presence of too many unknown 
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variables. Therefore, further simplifications are made. First, the phase of each 
individual signal (6.34) is assumed to be equal 


PA, q = PA V qd € [1,0] (6.36) 


which will later be ensured by a calibration prior to the actual measurement. 
Second, due to the horizontal and vertical alignment of the FEs, the angle 
variables can be partially equated to form a left, right, bottom and top angular 
position 


0 = 0) = 43 (6.37) 
Ox = 02 = 04 (6.38) 
Ya =Y41 =p (6.39) 
YT = 3 =W4 (6.40) 


Subsequently, the individual signal attenuations can be substituted by the prod- 
uct of two fictional attenuations that correspond to their angular positions (left, 
right, bottom and top) 


A, = Ap: AL (6.41) 
An = Ap: Ar (6.42) 
A3 = Aq: AL (6.43) 
Ag = Ar: AR (6.44) 


Additionally, for the sake of clarity, the following substitutes for the sinc- 
functions will be employed 


h¿(a) = sinc (= (sin(@,) - sin(a)) - dr (6.45) 


hg(B) = sinc (= (sin(%,) — sin(B)) - da] (6.46) 
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By inserting all these substitutes, (6.35) can be simplified to 


Xala, B) = NANN, exp {jpa} (6.47) 
- [ALAs - hL (a) - hg(8) + ArAg : hr(a) - hg (B) 
+ ALAT - hi (a) - hr(ß) + ArAr - hg(a) - hr(ß)] 


Factorizing the fictional attenuations and their corresponding sinc-functions 
allows to reformulate the equation as 


Xa (æ, B) = NANN, - exp {jpa} (6.48) 
- [AL - h(a) + Ar + hr(a)| - [Ag - he (B) + Ar: hr(ß)] 
In this form, the amplitude relations for the control of the azimuth and elevation 


AoA of the composed virtual target can be found individually by setting the 
partial derivative according to œ and £ equal to zero respectively 


Vxa(a, p) =0 a) 
Palen) = NANN, - exp {iva} (6.50) 
Vaz ande) i Aes Ohp(@) | [Ap - ha(ß) + Ar: hr(ß)] =0 
a da 
un = WENN; - exp Gea} (6.51) 
ol ol 
- [AL hy (a) + Ar: hg(a)] - | Ap: E +AT: > =0 


According to the zero-product property, a product is zero if one of its elements 
equals zero [DF03], which can be utilized to make the following extractions 


, Oh, (a) Perag Ohr (a) -0 


L R (6.52) 
ða da 
Ohg(fB) Ohr(B) _ 
Ag: OB + AT: OB =0 (6.53) 
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From this, the amplitude relations to control the horizontal and vertical position 
of the superimposed virtual target can be derived 


Ohg(a) 
AL _ da 
Tee (6.54) 
da 
Ohr(B) 
A ONT(B) 
Bu ee (6.55) 


Together with the adjusted amplitude-based superposition of the signals 


2 
Q 
Da Aq) _ (Ay + Ap + A3 + Ay)” 


o~P= 
Ro Ro 


(6.56) 
the required RTS channel attenuations for a given azimuth ag; and elevation 
Pig angle can be determined unambiguously. 


For reasons of completeness, the derivatives of the sinc-function substitutes are 
given as 


da mia - dyx + (sin(@,) — sin(a))? 
cos(a) cos (74 - (sin(d,) — sin(a)) - dix) 
a I o (6.57) 
sin(#,) — sin(a) 
Ahg(B) cos(ß) sin (x Mu (sin(4¿) — sin(ß)) - dix) 
OB Nt di (sin(Wq) = sin(ß))? 
cos(fB) cos (x a - (sin(w,) — sin(f)) - da) 
= (6.58) 


sin(W,) — sin(ß) 
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Measurement 


In order to verify the analytical expressions above, a measurement campaign 
was conducted using the same RTS system and configurations as in section 6.6. 
Four FEs were positioned at a distance of R¿ ~ 1m in front of the RuT in a 
semi-circular formation in the horizontal and a squared formation in the vertical 
plane. A metal sheet facilitated the placement of the FEs and shielded off most 
of the unwanted static radar reflections of the background. A photograph of the 
setup can be seen in Fig. 6.22. The angular positions of the FEs were measured 
individually with the RuT and are stated in Table 6.3. The RuT was again 
configured according to the parameters given in Table 6.1 with the exception 
that this time all Nix = 3 transmit antennas with a vertical element spacing of 
dix = 1/2 were used. 


Figure 6.22: Photograph of the two-dimensional arbitrary angle of arrival measurement setup O 
2022 IEEE [DNGZ22] 
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Table 6.3: RTS FE Angular Positions 


Frontend Azimuth Elevation 


1 -5.4° —8.8° 
2 4.5° =. 
3 —3.4° 8.4° 
4 3.8° 9.9 


As mentioned before, a calibration measurement was executed in order to 
achieve phase coherency among all RTS channels utilized. The process was 
similar to the one described in chapter 6.5, but adapted for the usage of four 
adjacent FEs. First, all channels were calibrated in terms of amplitude and 
range by performing single measurements with only one FE active at a time. 
Next, phase coherency was established sequentially in tuples with only two 
channels active simultaneously. For this, the angle error of the virtual target 
composed by the two channels was monitored while the delay of one of the 
channels was swept and the delay of the other was kept constant. The results 
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Figure 6.23: Measurement for the the phase calibration with the angle error over the artificial delay 
offset ATrts for different RTS FE tuples © 2022 IEEE [DNGZ22] 
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of these measurements are shown in Fig. 6.23, where the angle error is plotted 
as a function of the simulative delay offset At; between the respective FE 
tuples. Phase coherency is achieved at the point of minimum angular error 
whose corresponding delay offset is used to calibrate the respective channels 
for the subsequent AoA measurements. 


For the actual two-dimensional arbitrary angle of arrival measurement series, 
a total of 30 individual measurements comprising of six different azimuth and 
five different elevation angles were performed. The required RTS channel 
attenuations were determined according to (6.54) and (6.55). Fig. 6.24 displays 
the measured and set-point value of the AoA of the superimposed targets, as 
well as their deviation and the position of the FEs. It can be assumed that the 
angular errors are caused by the vertical and horizontal misalignment of the 
FEs and an imperfect phase calibration. Establishing phase coherency among 
four RTS channels simultaneously presents a challenging task as the calibration 
process can only be performed metrologically in pairs, which is likely to fall 
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Figure 6.24: Results of the 30 measurements showing the detected and set AoA, as well as their 
deviations and the angular position of the FEs © 2022 IEEE [DNGZ22] 
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Figure 6.25: Simulated and measured (a) azimuth and (b) elevation average angle error O 2022 
IEEE [DNGZ22] 


victim to propagating errors. For the calibration only vertically or horizontally 
directly adjacent pairs can be selected, since in this either the azimuth or the 
elevation angle error is evaluated. In addition, the multi-path reflections caused 
by the mechanical structure can also lead to deviations of the measured angle 
from its nominal value as has been shown before. 


Fig. 6.25 displays the average azimuth and elevation angle errors over the re- 
spective set-point angles. The simulated values serve as a reference of the 
expected error due to the inaccurate vertical and horizontal alignment of the 
FEs. This misalignment contradicts the simplification constraints made in the 
calculation of the channel attenuations that are required to create a composed 
target at the intended AoA. However, the simplifications are necessary to en- 
able an analytically unambiguous solution for the calculation of the channel 
attenuations, and failure to comply with their assumptions inevitably leads to 
angular errors. This source of error was also accounted for in the simulation. 


The measurement results verify the successful extension of the concept to 
superimpose the radar echoes of four vertically and horizontally distributed RTS 
channels in order to create virtual radar targets at arbitrary, two-dimensional 
AoA in both the azimuth and elevation domain simultaneously. Furthermore, 
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the results substantiate the necessity for an accurately executed positioning and 
calibration of the RTS FEs and the corresponding channels. 
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In recent years, the development of autonomous driving has evolved to new lev- 
els of sophistication. In order to guarantee proper operation and thus the safety 
of traffic participants, autonomous driving functions and the radar sensors on 
which they rely must be validated thoroughly. Since mere software simulations 
are not able to reliably cover all system and hardware characteristics, and road 
tests are highly labor-intensive and not reproducible, novel validation paradigms 
are required. 

For these reasons, radar target simulators (RTSs) are on the rise, as they pro- 
vide testing capabilities to validate radar sensors in-place and under laboratory 
conditions. They deceive a radar under test (RuT) by creating an artificial en- 
vironment comprising of virtual radar targets. For this, RTS systems receive, 
modify and re-transmit the radar signal emitted by the RuT. The simulated range 
and cross-range of a virtual target can be controlled by applying a delay to the 
signal and steering its re-transmission angle, respectively. 


To fulfill the task of reliably validating autonomous driving functions, RTS 
systems must be able to create credible, high-precision virtual radar targets 
at arbitrary positions. The accuracy with which the position of a target can 
be controlled was so far limited by the delay and angle rasterization, that is, 
the smallest delay increment and the angular spacing of the front ends (FEs) 
of the RTS. The concepts presented in this work lift these limitations and 
enable the generation of high-precision virtual radar targets at an arbitrary, 
continuous range and cross-range while simultaneously reducing the hardware 
effort and associated costs. These capabilities are needed in order to deceive 
modern imaging radar sensors that perceive real-world objects, such as cars 
and pedestrians, as point clouds of reflections, so that they cannot distinguish 
between real world and simulated virtual environment. 


In the classical approach, the cross-range and thus the angle of arrival (AoA) 
of a virtual target is simulated by multiplexing the re-transmitted radar signal 
between lateral distributed FEs, which only allows for targets at discrete angular 
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positions, inevitably leading to non-negligible deviations between set-point and 
actual AoA. In this work, a new technique to superimpose the radar echoes of 
adjacent FEs enables the generation of virtual radar targets at an arbitrary AoA 
that can be controlled by the amplitude attenuations of the corresponding RTS 
channels. It helps to reduce the number of FEs and channels required and 
thereby cuts the hardware cost, while at the same time drastically decreases 
the angular error and increases the accuracy with which the virtual targets are 
created. The presented concept overcomes the rasterization and discontinuity 
limitations and allows to generate artificial complex, multi-reflection traffic 
scenarios as needed for the validation of autonomous driving functions. It has 
also been expanded for the application of two-dimensional angle simulation 
being able to stimulate 4D imaging radars by creating targets that are distributed 
in both the azimuth and elevation domain simultaneously. 


Furthermore, this work enables the analysis of the performance of RTS systems 
that implement the above technique in regard of their angle simulation capa- 
bilities. It provides an estimation of the expected angle error as a function of 
the FE spacing and RTS channel imbalances in amplitude and phase. A signal 
model has been developed that grants an in-depth comprehension of the virtual 
target generation process of RTS systems. It offers guidance in the design cycle 
of high-precision, angle simulating RTS systems and supports the placement 
and distribution of the RTS FEs based on the desired angle accuracy and the 
performance parameters of the utilized RF components. In addition, the model 
allows to estimate the number of FEs and their spacing required given the RTS 
channel imbalances and the angular error tolerable. It could be shown that 
due to the coherent characteristic of the superposition of adjacent radar echoes, 
the spacing of the FEs may exceed the intuitively assumed limit of the angu- 
lar resolution of the RuT and still create a superimposed virtual target at an 
arbitrary AoA. Moreover, a calibration process has been developed that helps 
to reduce phase and amplitude mismatches among the RTS channels utilizing 
only a commercial radar sensor and its output. It circumvents the need for a 
raw data extraction which makes it suitable for the validation of sealed sensors 
in-place and as-is. 


The simulation of the range of a virtual target is controlled by the delay that 
is applied to the radar signal. In the classical approach, the smallest delay 
increment determines the rasterization at which the range of a target can be 
adjusted, and thus the range accuracy of the RTS. It must satisfy not only 
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the range resolution capability of modern radars, that are on the verge to be 
improved using greater bandwidths, but also the even finer range accuracy of the 
sensors. A new technique presented in this work enables a stepless adjustment 
of the delay so that virtual targets can be created at arbitrary range without 
the restrictions of a rasterization. This also bears the potential to simulate the 
movement of virtual targets within the measurement cycle of a RuT allowing 
to create non-ideal targets that feature realistic characteristics such as range 
and Doppler migration. This is achieved by using fractional delay filters which 
enable the continuous, stepless adjustment of the delay applied to the radar 
signal by controlling the filter coefficients. A signal model was developed 
that introduces a method to generate virtual targets that exhibit both range and 
Doppler migration characteristics by solely fine-tuning the simulated range at 
a specific update rate. 


In the course of this work, a digital RTS system has been designed and imple- 
mented that enabled the verification and analysis of the techniques mentioned 
above. It consists of modular front ends that operate in the automotive fre- 
quency band of 76 — 81 GHz and a digital back end based on an FPGA with 
integrated analog-to-digital (ADCs) and digital-to-analog converters (DACs) 
on which the superposition of adjacent radar echoes and the fractional delay 
filters have been realized. With this, the test setup allowed the generation 
of virtual radar targets at arbitrary lateral and longitudinal positions and with 
range Doppler migration characteristics. The obtained measurement results are 
in good agreement with those of the signal models, which verifies not only 
the applicability of the techniques, but also the estimation capabilities of the 
models. 


Both concepts presented in this work operate independently of the modulation 
scheme of the RuT employed, as they are both based on true time delay and 
phase shifting methods that mimic the signal propagation in the real world. 
They do not require access to raw data and can therefore be implemented for 
commercially available automotive radar sensors. Moreover, the concepts are 
not limited in regards of the methodology of the RTS, as they are applicable to 
both analog and digital systems. 
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Automotive radar sensors play a major role in the«ur- 
rent develöpment of autonomous driving, asthey are 
capable of detecting objects even under advefse weather 
conditions. Their functional operation must be validated 
in-place and under laboratory conditions, which is why 

a fully integrated test system is. required. Radar target 
simulators (RTSs) deceive a radaf.under test (RuT) by 
creating an artificial environment consisting of virtual 
radar targets. However, due’to the design of current 
RTS systems, virtual targets have so far been subject 
to rasterization in the range and angle of arrival (AoA) 
domain. In this work, new techniques are presented that 
overcome this deficiency and enable the generation of 
virtual targets at arbitrary high-precision positions.\This 
allows for continuous Movement of the targets and thus 

a more credible simulation environment. \ 
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is in the field of radar target simulation, in particular 
the high-precision generation of virtual radar targets. 
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